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PREFACE 


The  essence  of  engineering  is  the  ability  to  apply 
theory  in  a  practical  situation.  This  statement  seems 
harmless  enough  but  in  reality  can  be  the  proverbial  "bucket 
of  worms";  it  is  a  very  difficult  task.  At  times 
engineering  becomes  more  of  an  art  than  a  science.  This  is 
evident  when  we  consider  numbers  used  to  make  the  theory 
match  reality,  technically  referred  to  in  engineering  terms 
as  "efficiency  factors"  or  "proportionality  constants." 

This  thesis  project  highlighted  some  of  the  difficulties  of 
aeronautical  engineering.  It  started  innocently  enough  as 
an  attempt  to  determine  an  optimum  roll  axis  for  the 
air-to-air  tracking  task.  This  required  adequately 
describing  the  dynamics  of  the  attacking  aircraft,  target, 
and  the  relative  kinematics  between  the  two.  The  project 
ended  with  efforts  to  adequately  describe  the  attacking 
aircraft  dynamics,  an  effort  that  was  initially  supposed  to 
be  a  minor  portion  of  the  overall  project. 

1  owe  debts  of  gratitude  to  many  people  and  would  like 
to  acknowledge  a  few  of  them  here.  First  of  all  to  Captain 
Jim  Silverthorn,  who  provided  invaluable  assistance 
throughout  the  project,  upon  numerous  occasions  explaining 
the  finer  points  of  the  science  part  of  engineering  and 
providing  enlightenment  into  the  art  of  engineering.  To  Dr. 
Calico,  for  his  assistance  in  setting  up  the  project. 
Finally,  to  my  classmates  of  USAF  Test  Pilot  School  Class 
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82B,  especially  my  test  management  project  group  members. 
Major  Steve  heaps.  Captains  Steve  Pitotti,  Roger  "Flash" 
Keith,  and  Ed  Solski  (Canada),  and  Lt.  Samuel  “Shmulik” 
Nehama  (Israel).  They  put  up  with  my  unique  style  of 
leadership  and  were  always  ready  to  help. 
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ABSTRACT 


v_s»,A  nonlinear  simulation  and  a  linear  model  are  developed  for 
an  air-to-air  tracking  task  using  the  A-7D  DIGITAC  aircraft.  The 
nonlinear  simulation  includes  the  aircraft,  flight  control 


system,  and  target  dynamics,  and  the  relative  kinematics  between 


n 


the^aircraft  and  khtf' target.  The  linear  model  represents  the 
linearized  lateral-directional  aircraft  dynamics,  the  heading 
angle,  and  the  linearized  lateral  relative  kinematics.  The 
models  were  used  to  describe  the  effect  of  depressing  the 
aircraft  roll  axis  on  the  air-to-air  tracking  task  using  a  fixed 
gunsight.  — ^ 

A  flight  test  program  was  conducted  to  determine  the  effects 


of  depressing  the  aircraft  roll  axis  on  the  lateral-directional 

(Ms**' 

dynamics  and  on  the  air-to-air  tracking  task.  v7^0pen  loop\ tests  a^ 
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were  dene;  to  determine  the  effect  of  depressing  the  roll  axis  on 


the  lateral-directional  dynamics.  -Performed-  were;  (spiral,  Dutch 
(  roll,  and  roll  mode  tests,  and  lateral  frequency  sweeps^y^The 
~Tkt  most  noticeable  open  loop  effect-tS?  depressing  the  roll  axis^was 
the  decrease  in  roll  mode  time  constant.  Comparison  of  the 
flight  test  and  nonlinear  model  data  showed  good  correlation  for 
both  the  mechanical  flight  control  system  and  the  standard 
control  augmentation  system.  The  linear  air-to-air  tracking 
model  did  not  agree  with  flight  test  data  due  to  lack  of  a 
control  system  in  the  model. ^ Size  constraints  of  the  linear 


analysis  program  used,  prevented  incorporating  the  control  system 


system  dynamics  into  the  model. 


Closed  loop  tracking  tests  were  used  to  evaluate  the  effect 
of  a  depressed  roll  axis  on  the  air-to-air  tracking  performance. 
The  results  of  the  evaluation  showed  that  an  effective  roll  axis 
depression  angle  equal  to,  or  slightly  more  depressed  than  the 
fixed  gupsight,  were  favored  by  the  project  pilots  for  the 
air-to-air  tracking  task. 


I .  INTRODUCTION 


Substantial  work  has  been  done  in  the  past  with  digital 
flight  controls  (References  1,2).  The  majority  of  the  work  was 
on  optimizing  the  flight  control  laws  for  specific  tasks.  One 
task  for  which  the  flight  controls  were  optimized,  was  that  of 
air-to-air  tracking  with  a  fixed  gunsight.  The  air-to-air 
tracking  task  has  traditionally  been  the  most  difficult  task  in 
aerial  combat.  It  requires  precise  fuselage  pointing  without 
adversely  exciting  gunsight/bullet  stream  dynamics.  The 
air-to-air  tracking  task  consists  of  two  subtasks  which  are 
usually  accomplished  simultaneously.  One  is  the  longitudinal 
subtask  of  acquiring  the  proper  lead  on  the  target.  The  second 
is  the  lateral-directional  subtask  of  getting  into  the  target 
turning  plane.  This  thesis  considers  the  lateral-directional 
subtask  and  attempts  to  analytically  describe  it. 


BACKGROUND 

A  YA-7D  DIGITAC  aircraft  modified  with  a  Digital  Flight 
Control  System  (DFCS)  has  previously  been  used  to  investigate  the 
effects  of  flight  control  law  modifications  on  aircraft  handling 
and  task  performance  (References  1,2).  This  has  led  to  the 
development  and  refinement  of  two  multimodes.  Precision  Attitude 
(PA)  and  Flight  Path  (FP).  The  PA  multimode  is  primarily 
concerned  with  controlling  aircraft  rates  for  precise  fuselage 
pointing  tasks  such  as  air-to-air  tracking.  The  FP  multimode  is 


primarily  concerned  with  controlling  the  aircraft  velocity 
vector.  In  the  longitudinal  sense,  the  PA  mode  can  be  thought  of 
as  a  pitch  rate  command  system  and  the  FP  mode  as  a  normal 
acceleration  (g)  command  system.  These  two  multimodes  are  not 
completely  separate,  as  the  PA  mode  blends  into  the  FP  mode  for 
large  amplitude  aircraft  maneuvers.  The  PA  mode  also  contains  a 
submode  with  provisions  for  depressing  the  effective  aircraft 
roll  axis  by  using  the  rudder  to  create  a  yaw  rate  proportional 
to  the  roll  rate. 

Previous  DIGITAC  tests  showed  better  tracking  results  could 
be  achieved  by  depressing  the  effective  roll  axis  using  the  PA 
mode.  Initial  tests  were  designed  to  determine  the  "optimum" 
roll  axis  using  a  fixed  depressed  gunsight  (Reference  3).  This 
"optimum"  roll  axis  was  based  upon  pilot  opinion  and  tracking 
performance.  The  results  of  the  initial  tests  indicated  the 
optimum  roll  axis  was  between  131  to  174  milliradians  depression, 
for  a  two  "g"  target  at  1500  feet  range.  Dynamic  gunsights  were 
then  introduced  .  (Reference  4).  The  tracking  errors  achieved  with 
dynamic  gunsights  were  about  three  times  worse  than  with  the 
fixed  sight.  This  led  to  further  testing  with  a  dynamic  sight  to 
determine  if  a  different  "optimum"  roll  axis  existed  for  the 
dynamic  sight.  Recent  tests  included  the  variation  of  other 
tracking  variables,  such  as  range,  to  determine  their  effect  on 
the  location  of  an  "optimum"  depressed  roll  axis  (Reference  5). 
The  test  results  indicated  that  the  optimum  depression  angle  for 
the  roll  axis,  was  at  an  angle  greater  than  the  sight  depression 


angle.  However,  the  results  did  not  isolate  any  one  axis  within 
this  range  as  the  "optimum." 

OBJECTIVE 

The  objective  of  this  thesis  is  to  develop  an  analytical 
model  of  the  air-to-air  tracking  task  and  compare  it  to  flight 
test.  If  the  air-to-air  tracking  task  can  be  modelled  as  a 
linear  closed  loop  control  task,  then  conventional  linear  control 
analysis  tools  could  be  used  to  identify  promising  areas  for 
flight  test  and  increase  the  understanding  of  the  flight  test 
results  achieved.  There  is  a  large  knowledge  base  in  the 
understanding  of  classical  frequency  domain  design  techniques 
which  could  be  applied.  If  the  air-to-air  tracking  task  could  be 
accurately  modelled  and  represented  on  a  locus  and/or  bode  plots, 
then  various  configurations  could  be  examined  to  determine  their 
effects.  The  results  could  be  used  to  guide  flight  test  efforts. 

APPROACH 

This  thesis  develops  a  nonlinear  simulation  of  the 
air-to-air  tracking  task,  performs  a  linearized  analysis,  and 
then  compares  these  models  against  flight  test  data.  The  models 
were  developed  starting  with  few  assumptions.  This  permitted  an 
estimate  to  be  made  of  the  errors  as  various  simplifying 
assumptions  were  later  introduced,  and  thereby  provides  a  means 


of  determining  the  validity  of- the  model.  To  accomplish  a 
successful  modelling  effort  for  the  air-to-air  tracking  task/  the 
following  steps  were  identified. 

1.  Develop  a  nonlinear  simulation  of  the  tracking  task. 

This  included  six  degree  of  freedom  nonlinear  aircraft  dynamics, 
the  flight  control  system,  target  dynamics,  and  the  relative 
kinematics  between  the  target  and  the  aircraft.  This  step  was  to 
establish  the  baseline  against  which  all  further  assumptions  were 
to  be  compared.  The  nonlinear  simulation  incorporates  fewer 
assumptions  compared  to  any  of  the  linear  simulations.  This 
becomes  important  when  attempting  to  describe  a  complex  task 
which  may  contain  important  nonlinear  effects.  The  nonlinear 
simulation  would  have  the  highest  probability  of  accurately 
describing  the  actual  task  dyanmics. 

2.  Describe  the  tracking  task  with  linearized  equations 
using  state  variables.  This  will  allow  the  use  of  conventional 
linear  analysis  tools  (ie.  Root  Locus),  in  analyzing  the  task. 

3.  Validate  the  models.  This  means  comparing  the  aircraft 
and  flight  control  model  data  with  flight  test  data.  This  was  to 
be  acconqslished  with  open  loop  (no  "pilot  in  the  loop")  tests. 

4.  Evaluate  the  ability  of  the  models  to  describe  the 
air-to-air  tracking  task.  This  is  the  final  step,  which  combines 
the  target  dynamics  and  relative  kinematics  with  the  validated 
aircraft  and  flight  control  system.  This  was  to  be  accomplished 
with  closed  loop  ("pilot  in  the  loop")  tracking  tests. 

To  limit  the  scope  of  this  project  to  a  managable  level, 
only  the  effect  of  depressing  the  roll  axis  of  the  aircraft  on 
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J  the  air-to-air  tracking  task  was  considered.  The  modelling 

effort  will  consider  the  aircraft  at  one  flight  condition  with 
fixed  target  range  and  simple,  cooperative  target,  maneuvers. 

The  target  range  was  fixed  at  1500  feet  since  that  is 
approximately  the  "normal"  range  for  air-to-air  gunnery,  giving  a 
bullet  time  of  flight  of  less  than  a  second.  The  fixed  range  in 
combination  with  a  fixed  depressed  sight  and  a  constant  g  target, 
results  in  a  constant  depression  difference  between  the  aircraft 
roll  axis  and  the  sight.  Thus,  the  effect  of  a  change  in  roll 
axis  depression  on  the  tracking  task  can  be  evaluated  without  the 
complexity  of  additional  variables. 

SIGN  CONVENTION 

9 

The  following  sign  convention  is  used  throughout  this 
project. 

1.  The  positive  axes  and  forces  are  out  the  nose  for 
the  X  axis,  out  the  right  wing  for  the  Y  axis,  and  down 
for  the  Z  axis. 

2.  The  positive  moments  are  pitch  up,  roll  right,  and 
yaw  right. 

3.  The  positive  control  surface  movements  are  elevator 
trailing  edge  down,  right  aileron  up,  and  rudder  trailing 
edge  left. 

4.  The  positive  pilot  control  inputs  are  stick  back,  stick 
right,  and  rudder  left. 
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II.  DEVELOPMENT  OF  THE  NONLINEAR  EQUATIONS 

The  first  step  in  performing  a  system  performance  analysis 
is  to  adequately  describe  the  system.  This  chapter  presents  a 
description  of  the  aircraft,  control  system,  and  air-to-air 
tracking  task  and  a  description  of  the  nonlinear  differential 
equations  used  to  model  these  dynamics. 

AIRCRAFT  DESCRIPTION 

The  A-7D  aircraft  is  a  subsonic,  single-seat,  land  based 
light  attack  airplane  powered  by  a  single  Allison  TF-41-A-1 
turbofan  engine  rated  at  14,250  pounds  thrust  without 
afterburning.  The  swept  wing  is  mounted  high  on  the  fuselage 
with  a  slight  anhedral.  The  wing  control  surfaces  include  plain 
sealed  inset  ailerons,  spoilers,  leading  edge  flaps,  and  single 
slotted  trailing  edge  flaps.  The  tail  has  a  swept  vertical  fin, 
rudder,  and  an  all  moveable  horizontal  stabilizer. 

A- 7  DYNAMICS 

The  flight  condition  selected  for  this  study  was  a  cruise 
configuration  at  15,000  feet  pressure  altitude  and  0.6  Mach. 

This  is  a  flight  condition  for  which  aircraft  data  was  readily 
available.  Aircraft  data  for  this  configuration  is  tabulated  in 

Table  1.  Some  of  the  primary  assumptions  used  in  developing  the 
equations  of  motion  were: 


1.  The  X  and  Z  axes  are  in  the  plane  of  symmetry  and 
the  origin  of  the  axes  is  at  the  center  of  gravity  of 
the  aircraft. 

2.  The  mass  of  the  aircraft  is  constant. 

3.  The  aircraft  is  a  rigid  body. 

4.  The  flat  earth  approximation  is  appropriate. 

5.  The  flow  is  quasi-steady. 

The  A-7  aircraft  was  considered  as  a  three  dimensional  body 
with  six  degrees  of  freedom,  three  translational  and  three 
rotational.  A  set  of  nonlinear  coupled  ordinary  differential 
equations  were  used  to  describe  the  aircraft.  Reference  7  by 
Roskam  was  used  as  the  source  of  the  equations.  The  equations 
are  shown  in  tables  2,  3 ,  and  4.  These  equations  were  written 
for  the  aircraft  body  axes.  Straight  and  level  flight  was 
assumed  for  initial  conditions. 


TABLE  I 

A-7D  Cruise  Configuration  Data 

(Reference  6) 

Altitude 

h 

15,000  feet  PA 

Mach  Mo. 

M 

0.6 

Weight 

W 

25,338  lbs 

Center  of  Gravity 

CG 

28.7%  MAC 

Dynamic  Pressure 

5 

300  slugs/sq  ft 

Wing  Area 

S 

375  sq  ft 

Wing  Span 

b 

38.73  ft 

Mean  Aero  Chord  (MAC) 

c 

10.84  ft 

Moment  of  Inertia 

*xx 

15,365  slugs  sq  ft 

Moment  of  Inertia 

Iyy 

69,528  slugs  sq  ft 

Moment  of  Inertia 

Izz 

79,005  slugs  sq  ft 

Product  of  Inertia 

1xz 

-1,664  slugs  sq  ft 

•L* 


- 


& 


TABLE  II 


NONLINEAR  AIRCRAFT  EQUATIONS 


A/C  FORCE  EQUATI ONS 


U  -  VR  -  (W+Wq)Q  -  g  sin  0  +  qSCx/m  +  T/m 
V  -  (W+W0)P  -  UR  +  g  sin  cos  0  +  qSCy/m 
W  «  UQ  -  VP  +  g  cos  fS  COS  0  +  qSCz/ro 


A/C  MOMENT  EQUATIONS 


-  I1  P  Q  -  I2  Q  R  +  I3  ?3  S  b  Cx  +  14  q  S  b  Cn 

-  -H  p  *  -  h  tp2-R2)  +  *7  9  S  c  Cm 

-  ra  P  Q  -  I,  Q  R  +  I10  g  S  b  Cl  ♦  I1X  §  S  b  Cn 


A/C  KINEMATIC  EQUATIONS 


if>  m  P  +  Q  sin  i>  tan  0  +  R  cos  tan  0 

0  ■  Q  cos  -  R  sin 

Vm  0  ®in  ^  /  cos  0  +  R  cos  /  cos  0 

•C  «  W  /  (U2  +  V2  +  W2)1'2 
/3  “S’  V  /  (U2  +  V2)1/2 


TABLE  III 


INERTIAL  TERM  EQPATIONS 


Cl-Clyy-W1..”71  ‘  1-,IX*2/I**I«>  3 

t1**1/1**1**  +  <Iir*Iyy^Ixx’3  ^  [1_<Ixx  ^xx1!*33 

U-<l„2/Ixx1««,J 

^*xe^*xx*ex3  ^  (1-*XE  ^1xxIxx3 

^xxx”xzz^  ^  Xyy 
Xxz  ^  *yy 

°*5  ryy 

Clx*2/lxxIzz  ~  (Iyy-Ixx  1  X«)3  /  C1’Xx*  ^XxxXzz3 
Ixz^Izz  ^1+^Izz"Iyy  ^  Xxx*3  ^  ^1_Ixz  ^xx^z^ 
*xxz^XxxXzz3  ^  lWxz  ^XxxXzz3 

d/i„)  /  U-ixb2/ixxi«) 


TABLE  IV 


AERODYNAMIC  DERIVATIVES 


'  CLo  +  CL-*<  +  <T.«  *C/2U  +  CLq  qC/2“  +  CUe 

'  CDo  +  k  CL2 
'  sin«<  -  CD  cos  ®< 

C  C  pb/2U  +  C  rb/2U  +  C  fa  +  C „r„  <5 

y/f  yp  yr  yaa  y*r 

'  -CD  sin  «<  -  CL  cos  oc 

cl/f/1  *  clp  pb/2u  +  clr  rb/2°  +  clS*  A  +  clfr  * 

cn-*  +  C.i  *C/2U  +  Cmq  qc/2U  +  Cmfe  & 

cn4  4  ♦  Cnp  pb/2u  *  cnr  rb'2u  +  CnSL  +  cnfr  * 
-qSC  /m 

qSCy/m 

1  /  ARe  e  »  Oswald's  efficiency  factor  -  0.8 


0.0194 


CONTROL  SYSTEM  DYNAMICS 


The  A-7D  DIGITAC  DFCS  duplicates  the  standard  A-7  Control 
Augmentation  System  (CAS)  and  in  addition  provides  two  pilot 
selectable  multimodes,  Precision  Attitude  (PA)  and  Flight  Path 
(FP).  The  PA  mode  is  used  to  depress  the  effective  roll  axis  of 
the  aircraft.  Appendix  A  gives  a  more  detailed  description  of 
the  flight  control  system.  Thus,  the  pilot  has  a  choice  of  three 
flight  control  setups:  mechanical  only,  mechanical  plus  CAS,  and 
mechanical  plus  multimode. 

The  A-7D  CAS  was  modelled  with  differential  equations  that, 
when  coupled  to  the  A-7  dynamics  equations,  provide  the  control 
surface  inputs.  The  model  includes  both  the  mechanical  and  CAS 
systems  from  the  pilot  force  input  to  the  control  surface 
deflections.  The  control  system  block  diagrams  are  shown  in 
Figures  1  and  2  (Reference  8).  It  was  necessary  to  include  the 
mechanical  system  since  the  CAS  is  only  partial  authority.  The 
CAS  is  connected  to  the  mechanical  system  through  a  series  servo 
(mechanical  adder)  and  provides  control  input  shaping.  The 
dynamics  of  the  rate  gyros  and  series  servos  were  not  considered 
significant  and  were  not  modelled. 

The  elevator  axis  contains  both  pitch  rate  and  normal 
acceleration  feedback  through  the  DFCS  from  rate  and  acceleration 
sensors.  In  addition,  bobweights  within  the  mechanical  system 
also  sense  pitch  acceleration  and  normal  acceleration  which 
provides  a  mechanical  feedback  of  these  variables  to  the  elevator 
surface . 


FIGURE  1:  MECHANICAL  AND  CAS  -  PITCH  AXIS 


MECHANICAL  AND  CAS  ROLL  AND  YAW  AXES 


The  roll  axis  CAS  has  roll  rate  feedback  through  the  DFCS  as 
well  as  pilot  lateral  stick  inputs.  In  addition,  the  mechanical 
path  contains  a  feel  spring,  viscous  damper,  and  roll  isolation 
servo,  all  of  which  has  a  significant  effect  on  the  roll  dynamics 
of  the  aircraft. 

The  mechanical  rudder  system  is  a  simple  gain.  The  yaw  CAS 
has  side  acceleration  and  yaw  rate  feedbacks  and  an  aileron  to 
rudder  interconnect  scheduled  with  elevator  position. 

The  PA  rudder  system  (Appendix  A)  has  side  acceleration  and 
derived  sideslip  rate  feedbacks  and  an  aileron  to  rudder 
interconnect  scheduled  with  angle  of  attack.  One  of  the 
principal  features  of  the  PA  mode  is  that  various  effective  roll 
axes  can  be  obtained.  This  is  done  by  modifying  the  sideslip 
rate  feedback  as  shown  below.  Estimated  sideslip  rate  is  derived 
from  the  equations 

$  -  p«c  -  r  +  (g/U)  sin  ^  cos  G  [1] 

m 

For  a  normal  coordinated  turn,  /S  should  be  zero.  Feeding 
this  back  to  the  rudder  creates  a  yaw  rate  into  the  direction  of 
the  roll.  This  tends  to  minimize  the  amount  of  sideslip 
generated  during  rolling  maneuvers,  by  effectively  depressing  the 
roll  axis  to  be  aligned  with  the  aircraft's  velocity  vector. 

That  is,  the  yaw  rate  combined  with  roll  rate  results  in  an 
depressed  roll  axis  according  to  the  approximation: 

/°  -  r  /  p  [2] 

where  /O  is  the  roll  axis  depression  angle  and  both  r  and  p  are 
measured  in  the  body  axes.  Depression  of  the  roll  axis  beyond 


the  velocity  vector  is  achieved  by  replacing  the  angle  of  attack 
term  in  the  above  equation  (which  represents  the  depression  of 
the  velocity  vector),  by  a  fixed  angle  term  (Kp*) ,  representing 
the  desired  depression  of  the  roll  axis.  For  example,  Kpe*  *»  0.1 
would  tend  to  produce  an  effective  roll  axis  0.1  radians 
depressed  from  the  fuselage  reference  line. 

AIR-TO-AIR  TRACKING  TASK 

The  air-to-air  tracking  task  is  one  of  positioning  the 
attacking  aircraft  behind  the  target  and  in  the  target  turning 
plane  with  the  proper  lead  so  as  to  place  a  gunsight  on  the 
target  and  maintain  it  there.  The  gunsight  approximates  the 
location  of  a  bullet  one  time  of  flight  after  being  fired  from 
the  attacking  aircraft.  The  gunsights  range  from  a  sight  fixed 
with  respect  to  the  aircraft  to  an  undamped  sight  which 
represents  the  true  bullet  stream  dynamics. 

There  are  generally  two  phases  to  the  air-to-air  tracking 
task.  First  is  the  gross  acquisition  phase  where  the  attacking 
aircraft  is  maneuvered  without  reference  to  the  sight  into  the 
general  position  necessary  for  tracking.  At  this  point  the 
aircraft  would  be  close  to  the  target  turning  plane, 
approximately  1500  to  2500  feet  in  trail  near  to  the  appropriate 
lead  angle.  Then  the  second,  or  fine  tracking,  phase  starts 
where  the  sight  is  now  referenced  in  attempting  to  correct  both 
the  lateral  and  longitudinal  errors  between  the  sight  and  target 
Pitch  corrections  are  used  to  zero  longitudinal  errors  and  roll 
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corrections  along  with  some  rudder  are  used  for  lateral  errors. 
Discussions  with  pilots  experienced  in  air-to-air  tracking ,  led 
to  the  impression  that  to  correct  for  lateral  errors,  the  pilot 
would  first  match  his  bank  angle  to  the  target,  determine  the 
lateral  error,  then  adjust  the  bank  angle  to  fly  into  the  target 
turning  plane.  The  sight  was  not  precisely  used  whenever  the 
attacking  aircraft  bank  angle  did  not  match  the  target's  due  to 
pendulum  effect. 

Pendulum  effect  occurs  whenever  the  aircraft  is  rolled  and 
the  sight  axis  is  different  than  the  aircraft  roll  axis.  When 
the  roll  axis  is  above  the  sight  (the  normal  situation)  ,  the 
sight  will  appear  to  move  away  from  the  target  when  the  aircraft 
is  rolled  towards  the  target.  This  is  depicted  as  adverse 
pendulum  effect  in  Figure  3.  In  this  example,  the  pilot  must 
roll  right  in  order  to  translate  over  to  the  target's  turning 
plane.  In  so  doing,  however,  the  gunsight  will  appear  to  move  to 
the  left,  initially  increasing  the  perceived  lateral  error.  This 
gives  the  pilot  a  false  indication  of  the  magnitude  of  the 
lateral  error  and  does  not  provide  a  cue  as  to  when  to  roll  out. 
This  can  lead  to  sight  oscillation  about  the  target  as  the  pilot 
tries  to  guess  when  his  roll  axis  is  in  the  target  turning  plane. 
The  opposite  effect  is  achieved  when  the  roll  axis  is  below  the 
sight.  In  this  case  the  sight  will  appear  to  move  toward  the 
target  as  the  pilot  rolls  toward  the  target.  Once  the  sight  is 
on  the  target,  the  pilot  keeps  it  there  and,  in  theory  at  least, 
it  will  lead  him  into  the  target  turning  plane.  A  third 
situation  occurs  when  the  roll  axis  is  at  the  sight.  Here,  for  a 
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fixed  sight,  there  would  be  no  pendulum  effect  and  the  sight 
could  be  flown  to  wherever  it  needed  to  be.  In  practice,  a 
computing  sight  would  take  gravity  drop  into  account  which,  would 
cause  the  sight  to  be  off  the  aircraft  centerline.  This  would 
cause  a  pendulum  effect  where  an  attempt  to  correct  a  lateral 
error  by  rolling  would  create  an  indication  of  both  a 
longitudinal  and  lateral  error. 

This  pendulum  effect  caused  by  the  location  of  the  roll  axis 
is  what  led  to  the  present  investigations  into  the  possible 
existence  of  one  "optimum"  roll  axis.  In  this  thesis,  only  the 
lateral  fine  tracking  task  with  a  fixed  sight  was  examined.  The 
longitudinal  task  was  eliminated  to  simplify  the  analysis.  The 
fixed  sight  was  used  because  with  a  dynamic  gunsight  it  is 
difficult,  if  not  impossible,  to  distinguish  aircraft  dynamics 
from  gunsight  dynamics  during  analysis  and  flight  testing. 

Since  the  thesis  is  dealing  with  a  fixed  sight,  the  task  of 
analytically  representing  the  air-to-air  tracking  situation 
became  simply  one  of  describing  the  motion  of  the  target  with 
respect  to  the  attacking  aircraft. 

The  target  was  represented  as  a  point  mass  performing  a 
constant  speed,  constant  "g",  turn.  The  following  equations 
represent  the  target  performing  such  a  maneuver,  where  VT  is  the 
target  speed,  and  the  angles  <f/^  and  ©T  describe  the  target 
direction.  Target  bank  angle  (|6T)  was  included  for  reasons 
discussed  previously. 
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*T  “  0 

m  An  r  sin  j>r  /  Vr  cos  67- 
©r  *  (AnT  cos  J^r  “  9  cos  ®r)  /  vr 

*r  “  *V 

An^.“  g  /  cos  ^7.  (for  a  level  turn)  [3] 


The  relative  kinematics  between  the  target  and  attacking 
aircraft  were  described  by  referencing  both  the  target  and 
attacking  aircraft  to  an  earth  fixed  reference  frame,  then 
transforming  them  back  up  to  the  aircraft  body  reference  frame. 
Figure  4  illustrates  and  Appendix  C  shows  the  equation 
development. 


NONLINEAR  SIMULATION 


This  chapter  presented  the  development  of  the  aircraft 
equations  of  motion,  control  system  dynamics,  target  dynamics, 
and  relative  kinematics  between  the  target  and  the  aircraft. 

These  equations  and  the  understanding  of  the  air-to-air  tracking 
task  presented  in  this  chapter  forms  the  basis  for  the  air-to-air 
tracking  simulation. 

The  simulation  uses  a  fourth  order  Runge-Kutta  with  an 
integration  interval  of  0.05  seconds  to  solve  the  differential 
equations.  The  0.05  second  interval  was  necessary  due  to  the 
fast  rudder  actuator.  The  simulation  went  unstable  with  a  larger 
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FIGURE  4:  RELATIVE  KINEMATICS 


interval  (0.1  second),  while  no  difference  in  results  were 
obtained  with  a  smaller  interval  (0.025  seconds).  The  entire 
simulation  was  coded  in  the  Pascal  language  to  allow  it  to  be 
executed  on  a  microcomputer.  The  flow  of  events  for  the 
simulation  are  shown  in  Figure  5. 

In  the  simulation,  the  angle  of  attack,  pitch  attitude,  and 
elevator  deflection  represent  perturbations  away  from  nonzero 
trim  values,  although  the  trim  values  are  assumed  zero,  whereas 
the  other  variables  (excluding  U)  have  true  zero  initial 
conditions.  This  is  not  really  a  problem  for  longitudinal 
maneuvering,  but  for  lateral-directional  maneuvering  the  problem 
becomes  especially  evident  in  the  V  equation.  There  will  always 
be  an  initial  W  velocity  in  the  body  axes  if  the  aircraft  is  at 
an  angle  of  attack.  The  generation  of  a  roll  rate  will  therefore 
create  a  value  which  will  cancel  the  value  created  by  the  yaw 
rate  and  the  aircraft  will  effectively  roll  about  the  aircraft 
velocity  vector  (i.e.  normal  coordination).  Since  a  zero  W 
initial  velocity  was  assumed,  that  does  not  happen  and  a  false 
side  velocity  is  generated  which  causes  problems  with  the  flight 
control  system  feedbacks  which  attempted  to  cancel  the  side 

e 

velocity.  The  solution  was  to  include  an  initial  W  term  in  the  V 
•  * 

and  U  equations.  In  the  V  equation  it  has  no  effect  as  long  as 
there  is  no  roll  rate,  but  creates  a  nonzero  term  once  a  rolling 
motion  starts,  which  is  the  desired  effect.  The  effect  of  the 
initial  W  in  the  U  equation  is  similar  for  pitch  rate. 
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The  thrust  is  assumed  to  be  fixed  at  the  correct  value  for 
the  trim  condition.  The  correct  value  is  when  the  thrust 
balances  the  total  drag.  The  CLq  in  the  lift  equation  is  not  the 
zero  angle  of  attack  lift  coefficient,  but  is  instead  the  lift 
coefficient  at  the  trim  condition. 

The  inputs  were  the  initial  values  for  the  state  variables 
(only  U,  VT,  and  r^  were  required,  all  others  were  assumed  zero), 
pilot  lateral  stick  and  rudder  force  and  the  starting  and  ending 
times  for  the  simulation  (generally  zero  start  time).  The 
outputs  were  all  28  state  variables,  normal  acceleration,  target 
normal  acceleration,  lateral  sight  error,  longitudinal  sight 
error,  side  acceleration  angle  of  attack,  sideslip  angle,  each  as 
a  function  of  time. 

The  output  variables  were  stored  in  a  disk  file  while  the 
program  was  executing.  A  separate  program  was  used  to  print  out 
16  of  the  variables  to  a  line  printer.  The  output  file  could 
also  be  used  to  generate  a  plot  of  any  output  variable  versus 
time  through  a  separate  plot  program.  The  plot  was  done  on  a  dot 
matrix  printer.  Appendix  B  contains  a  listing  of  the  simulation 
program  and  print  program. 


III.  DEVELOPMENT  OF  LINEAR  ANALYSIS 


The  development  of  the  nonlinear  simulation  provides  the 
baseline  against  which  the  validity  of  a  linear  model  can  be 
checked.  The  linear  model  provides  the  proper  representation  of 
the  air-to-air  tracking  task  for  using  linear  control  analysis 
tools  (i.e.  Root  Locus).  The  linear  analysis  program  used 
(TOTAL,  Reference  9),  could  only  handle  up  to  a  tenth  order 
system  when  using  the  state  variable  setup,  thus  limiting  the 
size  of  the  linear  model.  The  nonlinear  model  was  28th  order. 
Therefore,  the  problem  must  be  simplified  more  than  by  the 
standard  linearizing  assumptions.  The  complexity  of  the  flight 
control  system  precluded  being  able  to  reduce  the  linear  system 
to  the  tenth  order.  Other  linear  analysis  methods  could  have 
been  used  which  would  have  allowed  a  higher  order  linear  system 
analysis,  however,  they  would  have  limited  the  ability  to  make 
changes  to  the  model  or  to  examine  other  outputs  from  the  model. 
The  flexibility  of  the  state  variable  approach  was  a  significant 
and  highly  desirable  advantage  although  the  size  constraints  of 
this  approach  proved  to  be  a  major  barrier  in  this  thesis. 

Although  the  aircraft  response  could  not  be  adequately 
represented  in  reduced  order  linear  form,  a  preliminary  linear 
analysis  of  the  air-to-air  tracking  task  portrayed  an  interesting 
phenomena.  The  following  development  of  the  preliminary  linear 
model  is  presented  as  an  introduction  to  the  phenomena  noted. 

The  preliminary  linear  model  represented  the  air-to-air 
tracking  task  in  state  variable  form.  Six  state  variables  were 


originally  used.  Four  were  the  usual  lateral-directional 
aircraft  state  variables  (v,  p,  r,  ^)  .  Additionally,  the  heading 
angle  (V0»  was  included  and  the  lateral  relative  kinematic 
equation  was  also  used. 

From  Equation  [4]: 

^2  =  +  ^r3  +  V|*  [(cos  5^r  cos  ©r)  (cos  f  sin  *5  sin  9  -  sin  V 

cos  *$)  +  (cos  ©rsin  fV)(sin  V  sin  sin  9  +  cos  ^  cos  £) 

-  sin  Oy  cos  9  sin  fi  ]  -  V 
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and  higher  order  perturbation  terms  can  be  neglected,  one 
obtains : 

r2  =  -  rrx  +  pr3  +  VT ( - p )  -  v  t 5D 

Combining  this  with  the  linear  equations  obtained  from 
Roskam  (Ref  7)  one  obtains  the  sixth  order  state  variable 
representation  of  the  lateral-directional  aircraft  dynamics  and 
lateral  relative  kinematics.  The  value  of  Nr '  in  these  equations 
was  artificially  increased  to  approximate  a  yaw  stability 
augmentation  system  (SAS) .  This  avoided  modelling  the  entire  yaw 
axis  control  system  with  all  its  associated  dynamics.  This 
produced  a  well  damped  (/=0.5)  dutch  roll  mode  as  would  be 
expected.  Next,  a  roll  attitude  feedback  loop  was  added  which 
approximated  the  inner  loop  dynamics  of  a  pilot,  since  most 
pilots  can  rapidly  close  a  roll  attitude  loop.  The  effect  of  the 


roll  attitude  loop  is  to  couple  the  roll  and  spiral  roots 
together  and  move  them  outside  of  the  dutch  roll  roots  (figure 
6).  This  results  in  the  preliminary  state  variable 
representation  of  the  air-to-air  tracking  task  which  now  takes 
the  form : 

X  =  [Ajx  +  [B]u 
m 

where  r  =  [v  p  r  ^  ^  rj]  and  u  =  Appendix  D  shows  the 

linear  model  development  and  the  A  and  B  matrices. 

This  linear  representation  does  not  include  any  control 
system  or  target  dynamics.  This  model  also  represents  the 
depressed  roll  axis  effect  by  positioning  the  target  above  and 
below  the  normal  aircraft  roll  axis  (adjust  the  r^  term).  A 
better  and  more  accurate  way  to  describe  the  effect  of  depressing 
the  roll  axis  on  the  air-to-air  tracking  task  is  to  include  the 
control  system  dynamics,  this  would  then  represent  the  actual 
means  of  depressing  the  roll  axis  and  also  the  effect  of  the 
control  dynamics  on  aircraft  handling  and  the  tracking  task.  But 
as  described  before,  limitations  in  the  control  system  analysis 
program  prohibited  this. 

The  results  of  the  preliminary  linear  analysis  are  shown  on 
a  root  locus  plot  in  Figures  6  and  7.  These  figures  represent 
the  closed  loop  poles  that  would  result  if  the  pilot  fed  back 
lateral  position  error  (^J  to  <fac.  Figure  6  represents  a  roll 
axis  elevated  50  milliradians  above  the  target  while  figure  7 
represents  a  25  milliradian  depressed  roll  axis  below  the  target. 
The  location  of  a  pair  of  complex  zeros  close  to  the  dutch  roll 
poles  effectively  eliminates  the  dutch  roll  in  the  lateral 


position  error  dynamics.  The  interesting  phenomena  noted  with 
the  effective  depression  of  thie  roll  axis  is  the  movement  of  a 
zero  on  the  real  axis.  It  starts  from  near  zero  on  the  positive 
axis  and  moves  to  the  right  as  the  roll  axis  is  depressed.  The 
zero  wraps  around  from  plus  infinity  to  minus  infinity  as  the 
roll  axis  reaches  the  sight  depression  and  moves  in  on  the 
negative  real  axis  as  the  roll  axis  is  depressed  further  than  the 
sight.  This  intuitively  appears  to  be  the  correct  movement  since 
a  zero  in  the  right  half  s  plane  indicates  an  initial  time 
response  in  the  opposite  direction  of  final  movement.  In  this 
case,  this  corresponds  to  an  initial  increase  in  lateral  error 
caused  by  the  adverse  pendulum  effect.  The  switch  from  positive 
to  negative  real  axis  also  occurs  when  the  sight  goes  from 
adverse  to  proverse  pendulum  effect. 

A  preliminary  hypothesis  was  formed  which  considered  the 
zero  location  as  an  indication  of  the  pilot  lead  required.  If 
the  pilot  were  to  react  as  a  simple  gain  in  feeding  lateral 
position  error  back  to  a  bank  angle  command,  then  the  72  mil 
depression  would  yield  a  lightly  damped  response,  while  the  147 
mil  depression  would  give  relatively  good  damping.  To  increase 
the  damping,  the  pilot  would  have  to  add  lead  to  draw  the  root 
locus  to  the  left.  The  further  out  on  t>»'  negative  real  axis  (or 
in  on  the  positive  axis)  the  roll  axis  depression  zero  was, 
indicated  more  pilot  lead  was  required.  Pilot  lead  was  equated 
to  pilot  compensation  required.  Increasing  the  compensation 
(increase  lead)  required,  increases  the  pilot  workload  and 
generally  results  in  poor  pilot  ratings. 


IV.  THE  FLIGHT  TEST  PROGRAM 


The  flight  test  program  was  designed  to  validate  the 
nonlinear  model  and  then  evaluate  the  ability  of  the  linear  model 
to  describe  the  air-to-air  tracking  task.  Five  depressed  roll 
axes,  equally  spaced  from  72  to  172  mils  depression  from  the 
fuselage  reference  line,  were  used  for  both  the  validation  and 
evaluation  phases.  A  fixed  gunsight  at  122  mils  depression  was 
used  for  the  tracking  evaluation.  Open  loop  (no  "pilot  in  the 
loop")  tests  were  done  to  validate  the  nonlinear  model.  Closed 
loop  ( "pilot  in  the  loop" )  tests  were  done  to  evaluate  the  linear 
model.  The  following  sections  contain  excerpts  from  the  flight 
test  report  (Reference  10)  and  explain  the  test  objectives, 
methods,  and  results. 


TEST  OBJECTIVES 

OPEN  LOOP  OBJECTIVES 

The  general  open  loop  testing  objective  was  to  provide  data 
to  determine  if  the  analytical  model  adequately  described  the 
actual  aircraft  lateral-directional  dynamics.  To  achieve  this, 
lateral-directional  characteristics  were  determined  for  the 
baseline  unaugmented  A-7D,  normally  augmented  A-7D,  and  each  of 
the  five  depressed  roll  axes.  Specific  tests  included 
determination  of  the  Dutch  roll,  roll,  and  spiral  mode 


characteristics.  In  addition,  the  frequency  response 
characteristics  of  the  DIGITAC  aircraft  were  evaluated  at 
specific  points  in  the  test  envelope  and  compared  with  linear 
model  derived  results.  The  results  of  the  open  loop  tests  were 
correlated  with  closed  loop  tracking  task  results  to  evaluate 
those  handling  qualities  related  factors  which  affect  the 
desirability  of  a  depressed  roll  axis  for  air-to-air  tracking. 

TEST  METHODS 

In  general,  standard  USAF  Test  Pilot  School  Plight  Test 
Techniques  (PTT's)  were  used  to  determine  the  lateral-directional 
characteristics  of  the  DIGITAC  aircraft  (Ref.  11).  All  tests 
were  initiated  from  a  trimmed  condition  at  15,000  feet  pressure 
altitude  and  350  KIAS.  These  FTT’s  are  summarized  below. 

1.  Dutch  Roll.  A  rudder  doublet  was  used  to  excite  the 
Dutch  roll  modes  of  the  aircraft  for  the  various  flight 
control  configurations.  In  addition,  aileron  doublets  were 
performed  to  determine  if  the  Dutch  roll  was  excited  by 
aileron  inputs. 

2.  Spiral  Mode.  The  aircraft  was  placed  in  a  stabilized  20 
degree  bank  turn.  At  test  initiation,  the  stick  was  released 
and  the  time  to  either  half  or  double  amplitude  measured. 


3.  Roll  Mode.  Aileron  rolls  were  performed  to  determine  the 
roll  mode  time  constant.  Step  aileron  inputs  were  used  to 
measure  roll  perfomance.  Roll  perfomance  was  measured  for 
both  left  and  right  rolls.  Full  deflection  aileron  rolls 
were  initiated  from  level  flight  with  the  input  held  through 
360  degrees  of  roll. 

4.  Frequency  Response.  Lateral  frequency  sweeps  were  used 
to  determine  the  aircraft  lateral  frequency  response. 
Sinusoidal  aileron  inputs  were  initat  sd  at  a  constant 
magnitude  and  frequency.  The  frequency  was  increased  in  steps 
to  the  maximum  frequency  obtainable  by  the  pilot. 

CLOSED  LOOP  OBJECTIVES 

The  closed  loop  testing  objective  was  to  evaluate  the  effect 
of  depressing  the  roll  axis  on  the  air-to-air  tracking  task.  To 
achieve  this,  each  of  the  five  depressed  roll  axes  were  evaluated 
several  times  by  each  pilot  with  a  pilot-in-the-loop  tracking 
task  using  a  maneuvering  T-38  airborne  target.  The  task  was 
designed  to  primarily  excite  the  lateral-directional  dynamics  of 
the  DIGITAC  aircraft.  The  performance  criteria  were  tracking 
within  three  mils  of  the  aimpoint  (intersection  of  T-38 
tailpipes)  for  four  consecutive  seconds  for  desired  performance 
and  one  second  for  adequate  performance.  This  criteria  was 
established  to  create  a  high  gain  tracking  task  which  should 


provide  better  differentiation  between  the  roll  axes. 

For  the  first  15  seconds  of  each  tracking  task,  the  T-38 
target  flew  a  level  3  g,  left  turn.  During  this  interval,  the 
DIGITAC  pilot  moved  the  pipper  from  an  initial  50  mil  lateral 
offset  and  attempted  to  precisely  track  (fine  track)  the  T-38. 
After  the  15  second  interval,  the  T-38  initiated  one  of  two 
maneuvers  to  cause  the  DIGITAC  pilot  to  further  excite  the 
lateral  axis  while  continuing  to  track  (maneuver  track) .  The  two 
maneuvers  were  a  loaded  S  turn  on  the  horizon  and  a  loaded  half 
barrel  roll  (reversal).  The  S  turn  provided  changes  in  target 
roll  direction.  The  reversal  provided  a  constantly  varying 
target  bank  angle.  The  pilot  was  not  told  which  depression  was 
being  evaluated  nor  which  of  the  two  maneuvers  the  target  would 
fly.  This  was  done  to  prevent  any  bias  in  the  evaluation  of  that 
roll  axis  depression. 

At  the  completion  of  the  maneuvering,  the  DIGITAC  pilot 
assessed  the  compensation,  workload  and  performance,  then  gave 
separate  Cooper-Harper  ratings  to  the  fine  tracking  and 
maneuvering  portions  (Ref.  12).  The  Cooper-Harper  rating  scale 
is  shown  in  Figure  8 .  Comments  were  made  and  recorded  both 
daring  and  immediately  after  each  manueuver.  Comments  especially 
elicited  were  those  concerning  predictibility  (pipper  position  vs 
control  input),  performance,  and  compensation  required  to  achieve 
tracking. 
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DEFINITIONS  FROM  TN-D-5153 


COMPENSATION 

The  measure  of  additional  pilot  effort 
end  attention  required  to  maintain  a 
given  level  of  performance  in  the  face  of 
deficient  vehicle  characteristics. 

HANDLING  QUALITIES 
Those  qualities  or  characteristics  of  an 
aircraft  that  govern  the  ease  and  preci¬ 
sion  with  which  a  pilot  is  able  to  perform 
the  tasks  required  in  support  of  an  air¬ 
craft  role. 

MISSION 

The  composite  of  pilot-vehicle  functions 
that  must  be  performed  to  fulfill  opera¬ 
tional  requirements.  May  be  specified  for 
a  role,  complete  flight  flight  phase,  or 
flight  subphase. 


PERFORMANCE 

The  precision  of  control  with  respect  to 
aircraft  movement  that  a  pilot  is  able  to 
achieve  in  performing  a  task.  (Pilot- 
vehicle  performance  is  a  measure  of 
handling  performance.  Pilot  perform¬ 
ance  is  a  measure  of  the  manner  or 
efficiency  with  which  a  pilot  moves  the 
principal  controls  in  performing  a  task.) 

ROLE 

The  function  or  purpose  that  defines  the 
primary  use  of  an  aircraft. 

TASK 

The  actual  work  assigned  a  pilot  to  be 
performed  in  completion  of  or  as  repre¬ 
sentative  of  a  designated  flight  segment. 


WORKLOAD 

The  integrated  physical  and  mental  effort  required 
to  perform  a  specified  piloting  task. 


FIGURE  8:  COOPER-HARPER  RATING  SCALE 


DATA  REDUCTION  AND  TEST  RESULTS 


OPEN  LOOP  RESULTS 

Data  were  reduced  from  telemetry  ground  station  time 
histories  of  the  test  aircraft  parameters.  Dutch  roll 
characteristics  were  determined  from  the  rudder  position  and 
angle  of  sideslip  traces.  The  spiral  characteristics  were 
determined  from  cockpit  instrumention  values  of  bank  angle  and 
time.  The  roll  mode  time  constants  were  derived  from  time 
histories  of  roll  rate.  The  frequency  response  of  the  DIGITAC 
aircraft  was  determined  from  time  histories  of  lateral  stick 
force  as  the  input  and  roll  rate  response  as  the  output  (Ref. 

13) . 

The  results  of  the  open  loop  tests  were  tabulated  for 
validation  of  the  analytical  aircraft  model.  The  dutch  roll  was 
heavily  damped  and  the  spiral  mode  nearly  neutrally  stable  for 
the  CAS  and  multimode  control  systems  as  expected.  In  addition, 
specific  test  parameters  were  compared  to  determine  any  trends 
which  might  occur  as  a  result  of  depressing  the  aircraft  roll 
axis.  The  parameters  selected  for  analysis  were  those  which 
would  impact  the  handling  qualities  and  performance  of  a 
depressed  roll  axis  aircraft  for  the  air-to-air  tracking  task. 
These  were: 

1.  Roll  Mode  Time  Constant.  Of  primary  concern  to  the  pilot 
during  the  air-to-air  tracking  task  is  the  time  required  for  the 


aircraft  to  reach  a  steady  state  roll  rate  in  response  to  a 
lateral  stick  input.  The  roll  performance  parameter  useful  in 
describing  the  roll  response  of  an  aircraft  is  the  roll  mode  time 
constant.  Physically,  it  is  the  time  for  the  aircraft  to  reach 
63%  of  its  steady-state  roll  rate  following  a  step  aileron  input. 
The  roll  mode  time  constant  directly  influences  pilot  opinion  of 
the  maneuvering  capability  of  an  aircraft. 

The  effect  of  depressing  the  roll  axis  on  the  roll  mode  time 
constant  is  shown  in  Figure  9 .  The  roll  mode  time  constant 
decreases  as  augmentation  is  added  to  the  basic  mechanical  flight 
control  system.  An  increase  in  roll  response,  evidenced  by  a 
decrease  in  ^  ,  is  observed  with  the  control  augmentation  system 
on.  This  is  due  in  part,  to  the  additional  aileron  input  from 
that  system  as  shown  in  Figure  10.  The  increase  in  roll  response 
is  also  due  to  a  rudder  input  from  a  CAS  aileron  to  rudder 
interconnect.  As  the  roll  axis  is  depressed  with  the  Multimode 
system,  the  roll  mode  time  constant  continues  to  decrease  as  a 
result  of  increasing  rudder  commands  used  to  effectively  depress 
the  roll  axis.  Pilot  performance  in  the  air-to-air  tracking  task 
should  increase  as  a  result  of  this  increased  roll  response. 

2.  Steady-State  Roll  Rate.  The  steady-state  roll  rate  of 
the  aircraft  was  shown  to  be  essentially  constant  for  all  roll 
axis  depressions.  This  result  was  expected  due  to  roll  rate 
limiting  factors  present  in  the  flight  control  system. 


FIGURE  10:  AILERON  DEFLECTION  COMPARISON  FOR  STEP  AILERON  INPUT 


3.  Steady-State  Yaw  Rate.  The  yaw  rate  occurring  during 
roll  maneuvers  increased  with  roll  axis  depression.  This  was  due 
to  larger  rudder  inputs  from  the  Multimode  system  increasing  roll 
axis  depression. 

4.  Lateral  Acceleration.  The  change  in  lateral 
acceleration  experienced  by  the  pilot  in  rolling  maneuvers  is 
shown  in  Figure  11.  This  figure  indicates  an  increase  in  lateral 
acceleration  for  rolling  maneuvers  with  the  Multimodes  engaged 
over  the  basic  aircraft  control  system.  The  lateral 
acceleration,  however,  appears  to  remain  relatively  constant  for 
the  first  four  roll  axis  depressions.  A  increase  in  lateral 
acceleration  is  noted  for  the  172  mil  depression  setting.  This 
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increase  may  impact  pilot  rating  as  a  function  pilot  comfort. 

The  difference  between  left  and  right  rolls  is  probably  ^ue  to 
the  instrumentation  pod  (asymmetric  store)  on  the  right  wing. 

CLOSED  LOOP  RESULTS 

Figures  12  through  14  depict  Cooper-Harper  ratings  given  by 
all  project  pilots  for  different  roll  axis  depressions.  These 
ratings  are  separated  for  the  fine  tracking  and  maneuvering 
tracking  tasks.  Each  pilot  flew  each  configuration  at  least  twice 
during  the  test. 

From  Figure  12  for  fine  tracking,  we  note  an  obvious 
tendency  for  better  ratings  toward  the  122  and  147  mil  roll  axis 
depressions.  Similiar  results  were  indicated  for  maneuvering 
tracking  as  shown  in  Figures  13  and  14. 

It  is  interesting  to  note  that  overall  for  the  fine 
tracking,  the  ratings  are  better  than  those  for  the  maneuvering 
tracking.  Pilots  found  it  was  more  difficult  to  track  the  T-38 
during  the  "S"  turn  and  reversal  than  through  the  level  3  g  turn 
fine  tracking. 

Pilot  comments  showed  definite  trends.  The  comments  clearly 
showed  pilot  dislike  for  the  72  mil  roll  axis  depression.  There 
was  not  a  single  favorable  comment  for  this  axis.  Although  the 
72  mil  roll  axis  depression  is  close  to  what  the  pilots  fly  with 
operations ’ ly,  they  found  that  the  more  depressed  roll  axes  did 
not  produce  as  much  sight  oscillation  and  therefore  made  lateral 
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COOPER-HARPER  RATINGS  BY  ALL  PILOTS 
FOR  MANEUVER  TRACKING  -  "S  TURN" 
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FIGURE  14:  COOPER-HARPER  RATINGS  BY  ALL  PILOTS 
FOR  MANEUVER  TRACKING  -  "REVERSAL" 


tracking  easier.  The  amount  of  favorable  comments  increased  as 
the  roll  axis  depression  increased  up  to  and  including  147  mils. 
Comments  for  the  172  mil  depression  were  all  distinctly  less 
favorable  than  for  the  147  mil  depression.  Pilots  felt  the 
lateral  response  was  too  quick,  which  tended  to  produce  lateral 
oscillations.  Selected  comments  typical  for  the  roll  axis 
depressions  are  listed  in  Table  5. 

There  were  some  notable  lessons  learned  during  these  pilot-in 
the-loop  evaluations.  One  was  that  the  performance  criteria 
should  reflect  only  changes  in  the  independent  variable:  project 
pilots  felt  that  some  of  the  Cooper-Harper  ratings  were 
artificially  lowered  by  errors  in  the  longitudinal  axis  when  only 
lateral  axis  changes  were  being  evaluated.  Another  was  that 
control  harmony  (pipper  movement  vs  control  deflection)  changed 
for  the  longitudinal  and  lateral  axes  for  some  roll  axis 
depressions.  This  caused  problems  with  pipper  predictability. 
Roll  axis  depression  preference  might  well  be  different  depending 
on  where  this  harmony  occurs. 
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TYPICAL  PILOT  COMMENTS  BY  ROLL  AXIS  DEPRESSION 


Anywhere  but  on  him 


72  MILS 


Not  predictable  in  performance 


Lateral  oscillations  with  aggressiveness 


Smooth  inputs  help 


97  MILS 


Borderline  predictability 


Some  lateral  oscillations 


Lateral  control  pretty  good 


122  MILS 


Almost  no  lateral  oscillations 


Never  got  off  of  fuselage  and  wings 


Good  predictability  -  very  stable 


147  MILS 


Easy  to  control  and  predict 


Good  predictability  and  responsiveness 


Some  tendency  to  overcontrol 


172  MILS 


Lateral  oscillations 


Lateral  predictability  is  a  problem 


TABLE  V:  TYPICAL  PILOT  COMMENTS  FOR  ALL  TRACKING  TASKS 
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V.  COMPARISON  OF  MODEL  AND  FLIGHT  TEST  RESULTS 


The  analytical  models  were  compared  to  flight  test  results 

for  aileron  rolls  and  frequency  sweeps.  The  nonlinear  model  was 

used  to  generate  time  histories  of  a  full  stick  deflection  step 

aileron  input.  Roll  rate  and  yaw  rate  were  the  primary  time 

histories  used  for  comparison  with  flight  test  time  histories. 

There  was  good  correlation  between  flight  test  and  model  results 

for  the  CAS  as  shown  in  Figure  15.  A  mechanical  system  only 

comparison  also  showed  good  agreement  (Figure  16).  This  level  of 

agreement  for  both  roll  and  yaw  rate  time  histories  indicate  that 

the  nonlinear  models  of  the  mechanical  and  CAS  control  systems, 

along  with  the  airframe,  closely  approximate  the  actual  aircraft. 

The  yaw  rate  match  could  probably  be  improved  by  modifying  the 

lateral-directional  cross  stability  derivatives,  primarily  Cyp  or 

C  t  . 
no  a 

A  check  of  the  roll  rate  to  lateral  stick  force  frequency 
response  against  a  frequency  response  plot  generated  from  the 
linearized  roll  CAS  model  shows  good  agreement  as  shown  in  Figure 
17.  The  roll  CAS  model  was  a  ninth  order  linearized  CAS  and 
mechanical  roll  flight  control  system  and  lateral-directional 
aircraft  dynamics  model  (Appendix  D).  It  did  not  include  the 
rudder  flight  control  system  nor  the  lateral  relative  kinematics. 
It  should  be  noted  that  the  flight  test  frequency  response 
characteristics  were  hand  reduced  hence  limiting  the  number  of 
discrete  points  and  thus  preventing  the  drawing  of  a  flight  test 
frequency  response  curve.  A  frequency  response  analysis  program 
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FIGURE  15:  MODEL  AND  FLIGHT  TEST  DATA  COMPARISON  OF  ROLL  AND  YAW 
RATES  FOR  STEP  AILERON  INPUT  -  MECHANICAL  AND  CAS 


YAW  RATE  -  DEG/SEC  ROLL  RATE  -  DEG/SEC 


TIME  -  SECONDS 

—  Flight  Test 
o  Nonlinear  Model 


FIGURE  16:  MODEL  AND  FLIGHT  TEST  DATA  COMPARISON  OF  ROLL  AND  YAW 
RATES  FOR  STEP  AILERON  INPUT  -  MECHANICAL  ONLY 


was  available  but  could  not  be  used  because  the  DIGITAC  data  was 
not  in  the  proper  format. 

There  were  two  differences  between  flight  test  and  model 
derived  results.  First  there  was  an  asymmetrically  loaded 
instrumentation  pod  on  the  right  wing  which  means  the  X  and  Z 
axes  are  not  exactly  in  the  plane  of  symmetry.  This  could 
account  for  some  of  the  differences  between  flight  test  and  model 
data.  The  second  difference  was  that  the  stability  derivatives 
used  in  the  nonlinear  simulation  and  linear  model  correspond  to 
Mach  =  0.60  while  the  flight  testing  was  done  at  Mach  =  0.68. 
However,  the  difference  in  dynamic  pressure  is  accounted  for 
since  stability  derivatives  were  used  in  the  nonlinear  model. 
Checks  of  stability  derivative  variation  with  the  difference  in 
Mach  numbers  (0.6  to  0.68)  showed  negligible  changes. 

The  results  achieved  during  closed  loop  testing,  generally 
agree  with  the  results  predicted  by  the  preliminary  linear 
analysis.  Pilot  comments  indicated  that  the  system  damping 
increased  as  the  roll  axis  depression  increased.  Pilot 
compensation  (pilot  lead  required),  as  shown  by  the  Cooper-Harper 
rating**,  decreased  as  the  roll  axis  depression  increased.  The 
results  obtained  at  the  most  depressed  axis,  172  mils,  differed 
significantly  from  what  was  expected.  The  expected  results  were 
a  further  decrease  in  pilot  compensation  required,  and  an 
increase  in  system  damping.  The  results  obtained  in  closed  loop 
tracking  tests  indicated  that  more  pilot  compensation  was 
required  and  that  the  system  damping  decreased  resulting  in 
lateral  oscillations.  This  indicates  the  possible  existence  of 


an  optimum  operating  point  or  region  on  the  root  locus.  This 
could  also  be  an  indication  that  other  factors,  such  as  the 
increase  in  lateral  acceleration,  are  having  a  significant  impact 
on  aircraft  lateral-directional  handling  qualities. 


VI .  CONCLUSIONS 


The  thesis  objective  of  developing  an  analytical  model  of 
the  air-to-air  tracking  task  and  comparing  it  to  flight  test  was 
accomplished.  A  nonlinear  simulation  and  a  linear  model  of  the 
air-to-air  tracking  task  were  developed.  A  comparison  of  the 
nonlinear  simulation  roll  and  yaw  time  histories  for  a  step 
aileron  input,  against  flight  test  data,  showed  good  agreement. 

A  linearized  roll  model  did  show  good  agreement  with  flight  test 
data  with  respect  to  frequency  response  for  roll  rate  to  lateral 
stick  force.  However,  model  size  limitations  of  the  linear 
analysis  program  prevented  the  accurate  linear  modelling  of  the 
air-to-air  tracking  task. 

The  results  achieved  during  open  loop  testing  characterized 
the  effect  of  depressing  the  roll  axis  on  the  lateral-directional 
dynamics.  The  most  noticeable  change  is  the  reduction  in  roll 
mode  time  constant  as  the  roll  axis  gets  further  depressed.  The 
closed  loop  handling  qualities  evaluation  determined  that  the  122 
and  147  mil  roll  axis  depressions  were  the  optimum  for  the 
air-to-air  tracking  tests.  These  results  are  valid  for  test 
conditions  of  350  KIAS,  15,000  feet  PA  with  a  3g  target  at  1,500 
feet  range  and  with  a  fixed  depressed  sight. 

RECOMMENDATIONS 

1.  Additional  work  should  be  done  on  the  nonlinear  model. 

The  multimode  control  system  should  be  incorporated  and  the 


tracking  loop  should  be  closed  with  the  addition  of  a  pilot 
model . 

2.  Further  linear  analysis  should  be  done  to  incorporate  the 
flight  control  system  effects.  The  investigation  into  the 
effects  of  depressing  the  aircraft  roll  axis  should  be  continued. 

3.  The  flight  test  data  reduction  procedure  for  the  DIGITAC 
aircraft  needs  to  be  modified  to  allow  the  use  of  computer 
routines  for  frequency  response  analysis. 
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DESCRIPTION  OF  FLIGHT  CONTROLS 


FLIGHT  CONTROLS 

The  standard  A-7D  flight  control  system  includes  a  control 
augmentation  system  (CAS)  implemented  by  feeding  stick  and  rudder 
pedal  forces  into  an  analog  computer.  This  CAS  is  in  addition  to 
a  mechanical  path  that  physically  links  the  primary  control 
surface  actuators  to  the  stick  in  each  axis.  Dual  electro- 
hydraulic  actuators  connected  in  series  with  the  primary  mech¬ 
anical  linkage  implement  pilot  inputs  with  additional  partial 
control  authority.  All  control  surfaces  are  fully  powered  by 
irreversible  dual  tandem  hydraulic  actuators.  Control  surfaces 
are  an  all  moving  horizontal  tail  for  longitudinal  control,  a 
conventional  rudder  for  directional  control,  and  a  combination  of 
ailerons  and  spoiler/deflectors  for  roll  control.  Analog  com¬ 
puters  in  the  standard  AFCS  were  removed  and  replaced  with  the 
multimode  (DFCS)  for  DIGITAC  testing  (3:179). 

MECHANICAL  FLIGHT  CONTROL  SYSTEM 


The  primary  mechanical  flight  control  system  remains 
unchanged  by  installation  of  the  DFCS.  Conventional  control  rods 
link  the  pilot's  center  stick  controller  to  control  surface 
actuators  in  both  longitudinal  and  lateral  axes.  Cables  link  the 
rudder  pedals  to  actuators  in  the  directional  axis.  Two 
independent  hydraulic  systems  power  the  control  surface  actuators 


during  normal  operation.  The  longitudinal  feel  system  consists 
of  a  dual  gradient  feel  spring,  two  bob  weights,  and  two  viscous 
dampers.  These  dual  bob  weights  produce  a  static  force  gradient 
proportional  to  normal  acceleration  in  addition  to  the  feel 
spring  force  gradient.  The  dual  bob  weights  also  produce  dynamic 
forces  proportional  to  pitch  acceleration.  With  viscous  damping, 
this  system  minimizes  the  possibility  of  pilot  induced 
oscillations.  The  longitudinal  trim  system  is  a  parallel  pulse 
type  using  an  actuator  in  parallel  with  a  dual  rate  feel  spring. 
Lateral  stick  feel  is  obtained  with  a  simple  linear  spring  and 
viscous  damper.  A  feel  isolation  actuator  in  the  lateral  linkage 
prevents  feedback  from  the  spoiler/deflector  breakout  device. 
Directional  feel  is  obtained  from  linear  springs  which  produce 
different  force  gradients  in  clean  and  landing  configurations. 
Directional  axis  trim  uses  the  standard  flight  control  system 
actuator. 

MULTIMODE  DIGITAL  FLIGHT  CONTROL  SYSTEM 

General: 

The  DFCS  is  a  dual  channel,  three-axis  system  using  whole  word 
computer  techniques.  It  provides  inner  loop  stability  and  con¬ 
trol  augmentation  functions  and  outer  loop  automatic  pilot  relief 
functions.  The  DFCS  was  designed  not  only  to  duplicate  the 
operation  and  functions  of  the  standard  A-7D  AFCS  (described  in 
chapter  II),  but  also  to  provide  pilot  selectable  control  modes 


(multimodes )  which  implement  optimized  control  laws  for  improved 
aircraft  operational  performance.  These  operational  modes 
essentially  provide  improved  control  of  either  aircraft 
flightpath  or  aircraft  attitude  depending  on  the  mode  selected. 
Multimodes  are  called  Pitch  and  Lateral  Flightpath  (FP),  and 
Pitch  and  Lateral  Precision  Attitude  (PA) .  Standard  A-7D  AFCS 
modes  are  selected  on  the  standard  AFCS  control 

panel  which  has  been  rewired,  while  multimodes  are  selected  from 
a  new  multimode  control  panel  in  combination  with  the  standard 
AFCS  control  panel.  The  DFCS  uses  existing  A-7D  AFCS  sensors  and 
actuators  in  addition  to  sensors  added  for  multimode  features  and 
for  dual  redundancy.  DFCS  hardware  consists  of  the  following: 

1.  Two  central  processor  units  (CPU's)  with  Honeywell 
HDC-301  processors. 

2.  Two  SEMS-8  memory  units  (8k  each). 

3.  A  multimode  control  and  display  panel  which  provides 
engage/disengage  switches  for  multimode 
functions,  engage/status  switches  for  CPU's,  and 
gain  change  capability  through  seven  three-position 


switches . 


4.  Two  remote  terminals  providing  interface  and  signal 
conditioning  for  both  aircraft  sensors  and  discrete 
switching.  These  terminals  also  provide  interface 
with  the  multiplex  data  bus  by  wire  or  fiber  optics. 

5.  Dual  rudder  linkage  position  transducers  and  a  second 
lateral  accelerometer. 

6.  Dual  aileron  position  transducers. 

CPU's,  memory  units,  and  interface  units  are  located  in  the 
left  avionics  bay  while  the  multimode  control  panel  is  mounted  in 
the  cockpit  on  the  left  outboard  console. 

In  addition  to  the  normal  A-7D  sensor  inputs,  the  DFCS  uses 
airspeed  and  altitude  inputs  from  the  air  data  computer,  angle  of 
attack  (AOA)  information  from  the  AOA  vane  on  the  flight  test 
nose  boom,  and  attitude  information  from  the  ASN-50  attitude 
reference  system.  The  noseboom  AOA  vane  is  the  primary  source  of 
AOA  for  multimode  control  law  computations.  For  safety  reasons, 
the  boom  vane  transducer  is  powered  from  the  DFCS  rather  than  the 
flight  test  instrumentation  power  supply.  Connection  of  the  boom 
AOA  transducer  to  flight  test  instrumentation  is  made  from  the 
DFCS  interface  unit  through  an  isolation  resistor.  This  config¬ 
uration  ensures  isolation  of  the  boom  AOA  signal  from  flight  test 
instrumentation  and  through  comparison  monitoring  provides  a  safe 
flight  control  input. 


Standard  A-7D  AFCS  Modes 


Control  modes  provided  by  the  DFCS  include  standard  A-7D 


AFCS  modes  of : 


1.  Yaw  stabilization 


2.  Pitch  and  roll  control  augmentation 

3.  Attitude  hold 


4.  Heading  hold 

5.  Heading  select 


6.  Altitude  hold 


Performance  and  function  of  each  of  these  modes  is  the  same  as 
described  in  the  A-7D  Flight  Manual  (14:1-91). 


DFCS  Optimized  Control  Laws: 


Stick- force-per-g  gradient  for  both  FP  and  PA  modes  is 
designed  to  match  the  midpoint  of  the  normal  A-7D  gradient.  In 
FP  and  PA  mode,  the  roll  and  pitch  stick  force  shaping  is  a 
parabolic  function  which  eliminats  the  dead  band,  produces  quick 
er  stick  command  response  to  initial  inputs  and  provides  smooth 
control  at  small  input  levels.  Gain  scheduling  varies  feedback 


gains  as  a  function  of  dynamic  pressure  in  both  roll  and  yaw 
axes.  The  yaw  axis  uses  kinematically  derived  sideslip  rate, 
measured  lateral  acceleration  and  pedal  position  as  inputs  in 


lieu  of  an  aileron-to-rudder  interconnect.  This  provides  im¬ 
proved  coordination  and  minimizes  sideslip  due  to  gusts.  The 
derived  signal  is  supplied  to  the  rudder  through  a  lag-lead 
shaping  network  to  enhance  directional  stability  at  high  AOA. 
Dynamic  pressure  scheduled  gain  compensates  for  varying  rudder 
effectiveness . 

Flightpath  Mode; 

The  FP  mode  provides  optimized  flightpath  response  for  large 
amplitude  combat  maneuvers  which  are  typical  in  air-to-ground 
bombing  or  air  combat  maneuvering.  Stick  force  sensors  are  used 
for  both  pitch  and  roll  inputs,  and  rudder  pedal  position  sensors 
are  used  for  yaw  inputs.  FP  mode  provides  rapid  normal  accel¬ 
eration  response  at  the  expense  of  increased 

pitch  rate  overshoot  and  rapid  roll  rate  response.  FP  mode  also 
includes  a  pseudo-neutral  speed  stability  which  decreases  stick 
force  required  to  maintain  trim  with  speed  changes. 

Precision  Attitude  Mode: 

The  PA  mode  provides  optimized  attitude  response  for  small 
amplitude  fuselage  pointing  tasks  which  are  typical  in  air-to-air 


and  air-to-ground  gunnery.  As  in  FP  mode,  stick  force  sensors 
are  used  for  both  pitch  and  roll  inputs  and  rudder  pedal  position 
sensors  are  used  for  yaw  inputs.  PA  mode  provides  optimized 
pitch  rate  response  for  small  perturbations  around  trim.  Gust 
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response  is  minimized  at  the  expense  of  normal  acceleration 
overshoots.  Pseudoneutral  speed  stability  is  implemented  as  in 
FP  mode.  The  PA  mode  in  pitch  has  an  automatic,  smoothly  blended 
reversion  to  FP  mode  type  characteristics  at  normal  acceleration 
perturbations  above  0.25  g's.  When  this  limit  is  exceeded,  the 
system  approximates  FP  response  characteristics  by  feeding  back 
normal  acceleration.  Roll  axis  control  laws  are  identical  to 
those  of  FP  mode.  The  yaw  axis  uses  kinematically  derived  side¬ 
slip  rate,  lateral  acceleration,  and  rudder  pedal  position  inputs 
as  in  the  FP  mode.  However,  in  PA  mode,  a  +0.05  g  dead  band  is 
inserted  in  the  lateral  acceleration  feedback  path  to  prevent  the 
system  from  opposing  initial  lateral  accelerations.  This 
provides  better  turbulence  response  as  yaw  rate  is  sensed  and 
damped,  resulting  in  lower  rates  with  less  weathercocking .  As  in 
the  FP  mode,  automatic  rudder  trim  is  provided,  and  dynamic 
pressure  scheduled  gains  are  used  in  both  roll  and  yaw  axes. 

An  additional  lateral-directional  submode  is  provided  in  PA 
mode.  This  second  submode  or  "gun"  mode  (PA  LAT  2)  compromises 
conventional  coordination  in  favor  of  a  pilot  selectable  depres¬ 
sed  roll  axis  as  long  as  lateral  accelerations  are  low.  It  makes 
the  aircraft  roll  about  the  depressed  roll  axis  rather  than  the 
flightpath  for  roll  rates  up  to  0.4  radians  per  second.  Should 
pilot  commanded  roll  rate  exceed  this  limit,  the  DFCS  will  limit 
control  surface  commanded  roll  rate  to  a  maximum  of  that  selected 
by  the  pilot  on  the  multimode  control  panel  or  0.4  radians  per 
second,  whichever  is  least. 
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Analytical  diagrams  of  the  multimode  pitch,  roll,  and  yaw 
axis  are  shown  in  Figures  A1 ,  A2,  and  A3  respectively.  Various 
configurations  are  selected  on  the  cockpit  control  panel  shown  in 
Figure  A4.  A  detailed  description  of  the  longitudinal,  yaw,  and 
roll  axis  is  presented  in  Reference  3. 
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FIGURE  Al:  PA  PITCH  AXIS 
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FIGURE  A3:  PA  YAW  AXIS 
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YAW  AXIS  DESCRIPTION 

Standard  A-7D  control  systems  are  complemented  in  the  yaw 
axis  with  a  stability  augmentation  system  (SAS).  Yaw-rate 
feedback  provides  damping  while  an  ARI  and  lateral  acceleration 
feedback  provide  turn  coordination  (1:260). 

In  contrast,  PA  and  FP  modes  provide  a  control  augmentation 
system  (CAS)  with  a  command  signal  from  rudder  pedal  position.  A 
combination  of  lateral  acceleration  and  yaw  rate  feedback  provide 
damping  and  turn  coordination  (Figure  A3).  An  ARI  scheduled  with 
AOA  rather  than  the  standard  horizontal  tail  position  was  intro¬ 
duced  for  DIGITAC  II  testing.  In  the  lateral  acceleration  loop, 
VjJ5  FP  and  PA  modes  are  identical  except  for  a  +1.6  ft/e<v?  deadband 

in  the  proportional  path  of  the  PA  mode.  This  path  allows  long 
term  lateral  trimming  in  both  PA  and  FP  modes,  and  improved 
lateral  damping  (weathercock  stability)  for  large  lateral  accel¬ 
erations.  Lateral  acceleration  response  differs  from  the  stan¬ 
dard  SAS  mode  in  that  lateral  acceleration  response  varies  with 
dynamic  pressure. 

LAT  1: 


To  provide  coordination  in  the  conventional  sense  (zero 
sideslip),  both  FP  and  PA  multimodes  operate  according  to  LAT  1 
submode  logic,  unless  LAT  2  is  selected.  The  LAT  1  submode 
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control  law  comes  directly  from  the  following  lateral  force 
equation : 


=  U0  (rs  +/2  )  -  g  sin  /  cos  0 


where:  a 

y 


U 

r 


o 

s 


/3 


/,  e 


lateral  acceleration  at  the  eg 

true  airspeed  along  longitudinal  axis 

stability  axis  yaw  rate 

sideslip  rate  relative  to  air  mass 

roll  and  pitch  attitude 


If  a^.  is  produced  mainly  by  sideslip,  maintaining  zero 
sideslip  rate  would  produce  zero  lateral  acceleration.  The 
criteria  for  /3  -  0  is  satisfied  by  keeping  control  error  small 

Since  yaw  rate  is  measured  in  the  body  axis,  stability  axis  yaw 
rate  is  obtained  by  resolving  yaw  rate  (r)  and  roll  rate  (p) 
through  AOA  ( oC ) ; 


r s  =  r  cos  c<  -  p  sin  <=< 


Assuming  a  small  angle  of  attack; 


And  adding  a  small  correction  term  because  the  yaw  gyro  has  a  2.5 
degree  tilt; 

rg  =  r  -  p«*  +  .  044p 

Now  assuming  a^  is  zero  because  is  small,  and  solving  for  ; 

/3  =  p  oc  -  r  -  .044p  +  (g  /  V)  sin f6  cos  0 

PA  and  FP  multimode  control  laws  implement  this  equation  as  shown 

* 

in  Figure  A3.  /3  is  the  control  element  and  is  kept  small  by 

driving  the  rudder  with  the  highest  loop  gain  possible.  The 
multimode  law  of  Figure  A3  schedules  Krm  with  lead-lag  shaping  to 
give  a  higher  gain  without  compromising  stability  margin  (7:33). 

LAT  2: 

An  additional  PA  submode  adds  a  roll  rate  to  rudder  feedback 

signal  to  the  A  equation.  LAT  2  control  laws  replace  the  p*  term 

in  the  A  equation  with  another  term,  K  (Figure  A3).  This 

term  drives  the  rudder  creating  a  new  "critical  yaw  rate"  which 

effectively  alters  the  aircraft  roll  axis  (7).  Figure  A5 

illustrates  the  vector  summation  of  body  roll  rate  and  roll  rate 

induced  body  yaw  rate  to  create  a  new,  depressed  roll  axis.  The 

amount  of  induced  yaw  rate  from  a  given  roll  rate  or  roll  rate  to 

rudder  gain,  is  a  function  of  K_ _  and  P3  setting  (Figure  A3).  A 

P 
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deliberate  stability  yaw  axis  uncoordination  is  allowed  in  favor 
of  altering  roll  axis  orientation. 


Note  that  there  is  a  limit  to  the  amount  of  uncoordination 
between  roll  axis  and  stability  axis.  When  the  limit  of  +  .4 
rad/sec  (25  deg/sec)  roll  rate  limit  is  reached,  no  further 
feed  forward  rudder  command  is  generated.  Additionally,  is 

no  longer  cancelled,  and  the  system  attempts  to  constrain  all 
further  sideslip  in  the  stability  axis. 

The  P4  switch  on  the  multimode  control  panel  (Figure  A4), 

sets  K  and  effectively  selects  reference  roll  axis.  Switch  P3 
poc  J 

applies  fixed  gains  of  1,  1.45,  and  2  to  the  reference  roll  axis 
which  can  further  depress  roll  axis. 

Table  A1 


Roll  Rate  To  Rudder  Gain 


Reference  Roll  Axis 


KpoC  Depression  WRT  Body  Axis 


0.072 

72 

mils 

0.097 

97 

mils 

0.122 

122 

mils 

0.147 

147 

milo 

0.172 


172  mils 


<** s+*) 

PROGRAM  A7; 

USES  TRANSCEND; 

CONST  Tdelta  =  0.05; 
Tint  =0.1; 


TYPE  VECTOR  =  ARRAY Cl.. 30 3  OF  REAL; 

VAR  XO, Xf inal ,F  :  VECTOR; 

Qbar , Qdot , adot , Antgt , Ptgt , Laz imuth, Lei evat i on, TO, Tf inal , Alpha, Beta, 
PILOTelev, PILOTai 1 , PILOTrud, Cx , Cl , Cm, Cn, An, Ay, Rdot , Pdot  :  REAL; 
ri,N, J, I, K, COUNT  :  INTEGER; 
b  :  CHAR; 

diskfile  :  STRING; 

OUTPUT  :  FILE  OF  RECORD 
STATE  : VECTOR; 

Normal accel ,Tgtaccel , azimuth, elevation, Sideaccel , 

AOA, Sidesl ip, Time  :  REAL; 

END; 


PROCEDURE  COEFFICIENTS (VAR  X:  VECTOR; adot  : REAL ;  VAR  Qbar, Alpha, 

Beta, Cx, Cl ,Cm,Cn, An, Ay  :  REAL); 

<*  This  procedure  computes  the  force  and  moment  coefficients  for  the  A7 
(*  DIGITAC  aircraft.  The  angular  units  are  radians  and  the  accelerat- 
<*  ions  are  in  ft/sec/sec.  The  sign  convention  is  positive  for  nose  up 
(*  nose  right,  and  roll  right.  The  positive  control  movements  are 
(*  elevator  trailing  edge  down,  right  aileron  up,  and  rudder  trailing 
(*  edge  left. 

CONST  k  =  0.0194;  (*  1/  pi  *  AR  t  Oswald's  efficiency  (80%)  *> 

mass  =  786.89; 

CIO  =  0.  174; 

Cl  alpha  =  4.412; 

Cl adot  =0.0;  (*  actually  Clalphadot  *) 

Cl q  =  1.0; 

Clde  =  0.596; 

CdO  =  0.0219; 

Cybeta  =  -0.7162; 

Cyp  =  0. 129; 

Cyr  =  0.3096; 

Cyda  =  -0.0502; 

Cydr  =  0.2006; 

Cl  beta  =  -0.0905; 

Clp  =  -0.346; 

Cl r  =  0.104; 

Cl  da  =  0. 1210; 

Cl  dr  =  0.0190; 

Cmalpha  =  -0.4636; 

Cmadot  =  -0.77;  ( $  actually  Cmalpbadot  *> 
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Cmq  =  -3.  "95; 

Cmde  =  -0.9056; 
Cnbeta  =  0.0722; 
Cnp  =  -0.00397; 
Cnr  =  -0.302; 
Cnda  =  0.0142; 
Cndr  =  -0.0917; 


chord  =■  10.84; 

(* 

feet  *) 

span  =  38.73; 

(* 

feet  *) 

Wingarea  =  375.0; 

(* 

square  feet 

rho  =  0.001497; 

(* 

15,000  ft.  1 

Clift, Cdrag, Cy,Cz  : 

REAL; 

*) 

MSL 


slug/cubic  ft.  *) 


(* 

FCi3  is  the  time  derivative 

of  X  C i 3 

*) 

(* 

XC  1 3  = 

U 

XC  173 

ss 

demanual 

*) 

(* 

XC23  = 

V 

XC  183 

= 

deelectric 

*) 

(* 

XC33  = 

W 

XC  193 

= 

de 

*) 

(* 

XC43  = 

P 

XC203 

= 

dalmanual 

*) 

(* 

XC53  = 

Q 

XC213 

= 

da2manual 

*) 

-  -  - 

(* 

XC63  = 

R 

XC223 

= 

da3el ectri c 

*) 

<* 

XC73  * 

Phi 

XE233 

=: 

da4electric 

*> 

(* 

XC83  = 

Theta 

XC243 

= 

da 

*) 

(* 

XC93  = 

Psi 

XC253 

= 

rdot  feedback 

*) 

(* 

XC 103= 

Target  Velocity 

XC263 

ss 

Ayl 

*> 

(* 

XC 1 1 3= 

Target  Psi 

XC273 

= 

Ay2 

*) 

(* 

XC 123= 

Target  Theta 

XC283 

= 

dr 

*> 

(* 

X  C 133  = 

Target  Phi 

*) 

(* 

XC 143= 

Relative  Pos.  - 

X 

*) 

(* 

XC 153= 

Relative  Pos.  - 

Y 

*) 

(* 

XC 163= 

Relative  Pos.  - 

Z 

*> 

vTv 

•A' 


BEGIN 

Qb ar :  =  (0.  5) *rho* (XC 1 3*XC 1 3+XC23*XC23+XC33 *XC33 ) ; 

Alpha: =XC33/SQRT (XC 1 3*XC 1 3+XC23 *XC23+XC33*XC33 >  ; 

Beta: =XC23 /SORT (XC 1 3 *X C 1 3+XC2 3 *X C23 >  ; 

Clift: =C1 0  +  Cl  alpha* Alpha  +  Cl adot * ( (adot Ichord) / (2*XC 1 3 ) >  +  Clq*( 

( X  C53  *chord) / <2*XC 1 3 ) >  +  Clde*XC193; 

Cdrag : =CdO  +  k*Cl i f t *C1 i f t ; 

Cx:  =  -Cdrag*cos (Alpha)  +  Clift*sin (Alpha) ; 

Cy:=Cybeta*Beta  Cyp*  (  (XC43 *span )  /  (2*X C  1  3  )  >  +  Cyr * ( (XC63*span) / (2*XC 1 3 > > 

+  Cyda*XC243  +  Cydr*XC283; 

Cz :  =  -Cdrag*sin (Alpha)  -  Cl i f t tcos (Alpha) ; 

Cl :=Clbeta*Beta  +  Clp* ( (XC43*span) / (2*XC 1 3 > )  +  Clr* < (XC63*span> / <2*XC 1 3) ) 
+  Clda*XC243  +  Cldr*XC283; 

Cm:  =Cmal  pha*Al  pha  Cmadot * ( (adot *chord ) / < 2* X r 1 3 > )  +  Cmq*  (  (  X  C53  ichord ) 
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"•j  -r .  r,  r,  7.  \r.  . . i . ■  .  >  v- 


/<2*XC13))  Cmde*XC193; 

Cn:=Cnbeta*Beta  +  Cnp* < <XC43*span ) / <2*XC 1 3 > >  +  Cnr* ( <XC63*span) / <2*XC 1 3) ) 
+  Cnda*XC243  +  Cndr*XC283; 

An:  =  -  (QbartWi  ngarealCz  )  /mass; 

Ay:  =  (Qbar *Wingarea*Cy) /mass; 

END; 


PROCEDURE  A7DYN (VAR  X:  VECTOR; Qbar. A1 pha, Cx , Cl , Cm, Cn, An, Ay  :  REAL; 

VAR  F:  VECTOR; VAR  adot  :  REAL); 

CONST  gravity  =  32.2  ; 

11  =  -0.03413; 

12  =  0.62048; 

13  =  0.000065232; 

14  =  -0.000001374; 

15  =  -0.9153146; 

16  =  -0.023933; 

17  =  0.0000143; 

18  -  -0.684845; 

19  =  -0.0341306; 

110=  -0.000001374; 

111=  0.00001269; 
mass  =  786.89; 
chord  =  10.84;' 
span  =  38.73; 

Wingarea  =  375.0; 

Thrust  =  3272.0; 

WO  =  40.0;  <*  Initial  AOA  -  trim  W  velocity  *) 


BEGIN 

<*  A/C  FORCE  EQUATIONS  *> 

FC13:=XC23*XC63  -  (XC33+W0) *XC53  -  gravi ty*si n ( X C83 )  + 
(Qbar*Wingarea*Cx > / (mass)  +  Thrust/mass; 

FC23 : = (XC33+WO) *XC43  -  XC13*XC63  +  gravi ty*si n (XC73 ) *cos (XCB3 )  +  Ay; 
FC33 :  =XC 1 3*XC53  -  XC23*XC43  +  gravi tyicos < XC73 ) *cos (XC83 )  -  An; 
adot;=FC33/XC13  -  Alpha* (FC 1 3/XC 1 3 ) ; 

(*  A/C  MOMENT  EQUATIONS  *) 

FC43 : =1 1 *XC43 *XC53  -  I2*XI53*XC63  +  I3*Qbar*Wi ngarea*span*Ci  +  14* 
Qbar*Wi ngarea*span*Cn ; 

FC53:=(-I5)*XC43*XC63  -  16* < X C4 3  *  X C4 3-X [63 *X C63 )  +  Qbar *Wi ngarealchord 
*Cm* 17; 

FC63 ;  =I8*XC43*XC53  —  I9*XC53*XC63  •*-  1 10*Qbar  *Wi  ngarea*span*Cl  +  Ill* 
Qbar*Wingarea*span*Cn; 

(*  A/C  KINEMATIC  EQUATIONS  *) 


FC73:=XC43  X C53 *si n < XC73 > *si n < XC8 3 ) /cos ( X C83 >  +  X[63*cos (XC73 > * 


'  .....  - .  ■  . 

%  • 


M2*XC13>>  +  Cmde*  X  C  1 9  3  ; 

Cn: =Cnbeta*Beta  +  Cnp * < < X C 4 3  *  span ) / < 2 *  X C 1 3 > )  4  Cnr * ( ( x C63 *span ) / < 2* X C 1 3 ) 
+  Cnda*XC243  +  Cndr*XC2S3; 

An:=  -(QbarlWingareaJCz ) /mass; 

A y:=  (Qbar *Wingarea*Cy > /mass; 

END; 


PROCEDURE  A7DYN (VAR  X:  VECTOR; Qbar . A1 pha , Cx , Cl , Cm, Cn , An , Ay  :  REAL; 

VAR  F:  VECTOR;  VAR  adot  :  REAL); 

CONST  gravity  =  32.2  ;  [ 

11  =  -0.03413; 

12  =  0.62048; 

13  =  0.000065232; 

14  =  -0.000001374; 

15  —  —0.9153 1 46; 

16  =  -0. 023933; 

17  =  0.0000143; 

18  =  -0.684845; 

19  =  -0.0341306; 

110=  -0.000001374; 

111=  0.00001269;  ! 

mass  =  786.89; 
chord  =  10.84; 
span  =  38.73; 

Wingarea  =  375.0; 

Thrust  =  3272.0; 

WO  =  40.0;  (*  Initial  AO A  —  trim  W  velocity  *) 


BEG  !  N 

(*  A/C  FORCE  EQUATIONS  *> 

PC  1 3:=XC23*XC63  -  (XC33+W0)  *XC53  -  gr avi ty * s i n < X C8 3 >  ■*- 
(QbartWi ngarealCx  )  /  (mass )  Thrust/mass; 

PC23:=(XC33  -*-WO )  *  X  C  4  3  -  XC13*XC63  +  t.r  avi  ty  *si  n  (  X C 73  >  *cos  ( X  C83  )  +■  Ay; 
FC33;=XC13*XC53  -  XC23*XC43  +  gr avi  t v *cos ( X C 7 3 ) *cos ( X C83 )  -  An; 
adot; =FC33/XC 1 3  -  AI pha* (FC 1 3 / X C 1 3 > ; 

<*  A/C  MOMENT  EQUATIONS  *> 

FC43  :  =1  1  *XC43  *XC53  —  I2*XC53*XC63  *■  I3*Qbar  tWi  ngareaispaniCl  ■*-  14* 
Qbar*Wingarea*span*Cn; 

P C 5 3  :  =■■  (  - 1 5)  *  X  C  4  3  *  X  C 6  3  —  16*  ( X  C  4  3  *  X  C  4  3  —  X  C63  *  X  C  6  3  )  ♦  QbarlWingareatchord 
*Cm* I 7 ; 

F 1 6  3  ;  =  1 8  *  X  C  4  3  *  X  C  5  3  —  1 9*  X  C  5  3  *  X  t  3  4  I  10*Qbar  *Wi  ngarea*span*Cl  Ill* 

Qbar  *Wi ng area* span  *Cn ; 

<*  A/C  KINEMATIC  EQUATIONS  * ' 


PC73:=XC43  *  X  C53  *si  n  (  X  t  7  3  )  *si  n  (  X  f  8  3  wCOs  (  X  [  8  J  »  4  X [ 6 3 *cos ( X C 7 3 > * 


v' 


si n ( XC83 ) /cos ( X  C83 ) ; 

FC83 : =X C53  *cos ( XC73 )  -  XC6 3 *si n ( X C7 3 ) ; 

FC93:=XC53*sin(XC73)/cos(XC83>  +  X C63 *cos ( X C7 3 ) /cos ( X C83 ) ; 
END; 

PROCEDURE  ASSIGN (VAR  F:  VECTOR;  VAR  Pdot,  Qdot,  Rdot  : REAL) ; 
BEGIN 

Pdot:=*FC43; 

Qdot : =FC53 ; 

Rdot: =FC63 ; 

END; 


PROCEDURE  TGTDYN (VAR  X:  VECTOR;  VAR  F:  VECTOR; 

VAR  Lazimuth, Lei evat i on, Antgt  :  REAL); 


CONST  gravity  =  32.2; 


VAR  Phi , Ph2, T1 , T2, Psl , Ps2, PsTl , PsT2, TT 1 , TT2, PhTl , PhT2  : 


REAL; 


BEGIN 

(*  TARGET  DYNAMICS 

Phi: =cos ( X  C  7  3 ) ; 

Ph2: =sin (XC73 ) ; 

T1 : =cos (XC83 ) ; 

T2: =sin (XC83 ) ; 

Psl: =co5 ( X  C  9  3 ) ; 

Ps2: =sin (XC93 ) ; 

PsT 1 : =cos ( X  C 1 1 3 ) ; 
PsT2: =si n(XC113) ; 

TT 1 :  =cos ( X  C 1 2  3 ) ; 
TT2: =sin ( X  C 123 ) ; 
PhTl : =cos (XC 133 ) ; 
PhT2:=sin (XC133) ; 


Antgt: =gravi ty/PhTl ;  (#  Assumin 

FC 103 : =0.0; 

Ft  1 1 3 :  = (Antgt *PhT2> /  (  X  t 103 *TT 1 > ; 

FC 123 := (Antgt *PhTl-gravi tytTTl > /XC103; 
FC 133: =Ptgt ; 


(#  Assuming  a  level  turn  *> 


RELATIVE  KINEMATICS 


fl 


FC143:=XC63*XC153  -  XC53*XC163  +  X C 103  * ( (TT1 *PsTl *T1 *Ps 1 )  +  <TTl*PsT2 
*Tl*Ps2)  +  (TT2*T2) >  -  XC13; 

FC153:=-XC63*XC143  +  XC43*XC163  ♦  X C 103 *(< (PsTl *TT 1 >* (Psl *Ph2*T2  -  Ps2 

*Phl>)  +  (  (TTl*PsT2> * (Ps2*Ph2*T2  +  Psl*Phl)>  -  (TT2*T1  *F'h2>  >  -  XC2 
FC163:=XC53*XC143  -  XC43*XC153  +  X C 103  * ( < (PsT 1 *TT 1 > * (Psl *Ph 1 *T2  +  Ps2* 
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Ph2) )  +  ( (TTl*PsT2> * (Ps2*Phl*T2  -  Ph2*Psl>)  -  ( TT2*T 1 *Ph 1 ) ) 
Laz i muth : =atan (XC153/XC143) ; 

Lelevation:=atan (-XC163/XC143) ; 

END; 


PROCEDURE  ELEVATOR (VAR  X:  VECTOR; An, Qdot „ PI LOTel ev  :  REAL; 

VAR  F;  VECTOR); 

CONST  gvslb  =  0.3; 

VAR  fc-fs,  del ,  decommand  :REAL; 

BEGIN 

(*  ELEVATOR  -  XC193  is  the  elevator  movement  (ie.  delta  e  >  *) 

IF  ABS (PI LOTel ev)  <  8  THEN 
kf s ;=ABS( PI LOTel ev) /32 
ELSE 

k-f  s:  =0.  25  +  (ABS  (PILOTel  ev)  -  8)  *0.344; 
del ; =PILOTelev  -  O. 00149* (An-32. 2)  -  0.167*Qdot; 

FC 173 : = (-0. 5*kf s) *XC 173  +  0.5*del; 

FC 181 : = (-XC 183  -  (An-32.2)  +  10*PIL0Telev) *1 . 818; 
decommand: =XC 173  +  (XC 183-460*XC53 ) *0. 00054; 

IF  decommand  >  0.262  THEN  decommand: =  0.262 
ELSE 

IF  decommand  <  -0.262  THEN  decommand: =-0. 262; 

FC 193 : =— 20*XC 193  -  20*decommand; 

END; 


PROCEDURE  AILERON (VAR  X:  VECTOR;  PILOTail  : REAL;  VAR  F;  VECTOR); 
VAR  dacommand  ;  REAL; 

BEGIN 

(*  AILERON  -  XC243  is  the  aileron  movement  (ie.  delta  a  )  *) 

FC203 : =—12. 8*XC203  +  2*PIL0Tail; 

FC21 3; =-12.5*XC21 3  +  1.24*XC203; 

FC223 : =-3. 33*XC223  +  O. 70*PIL0Tai 1 ; 

FC233 : =-10*XC233  +  10*XC223; 
dacommand:=  (XC233-XC43) *0. 1 ; 

IF  dacommand  >  0.175  THEN  dacommand: =0. 175 
ELSE 

IF  dacommand  <  -0.175  THEN  dacommand:=  -0.175; 

FC24  3;  =— 20*XC243  +  20*dacommand  +  20*XC213; 


PROCEDURE  RUDDER (VAR  X:  VECTOR;  PILOTrud,  Ay,  Rdot  :  REAL; 

VAR  F:  VECTOR) ; 

VAR  kl,  drl,  dr 2,  drcommand,  drpilot  :REAL; 

BEGIN 

(*  RUDDER  -  XC283  is  the  rudder  movement  (ie.  delta  r  )  *> 

drl:=0-0;  (*  Initialize  drl  *> 

(*  14  lb  breakout  force  t) 

IF  ABS (PILOTrud)  <  14  THEN  PILOTrud: =0. 0 

ELSE  IF  PILOTrud  >  14  THEN  PILOTrud : =PIL0Trud-14. 0 

ELSE  IF  PILOTrud  <  -14  THEN  PILOTrud: =PILOTrud  +  14.0 


drpilot: =0. 001064*PIL0Trud; 

it  6  degree  pilot  rudder  command  limit  in  cruise  *) 

IF  drpilot  >  0.105  THEN  drpi 1 ot : =0. 105 

ELSE  IF  drpilot  <  -0.105  THEN  drpi 1 ot : =-0. 105; 

(*  ARI  is  scheduled  with  elevator  (XC193).  *> 

(*  Normally,  the  ARI  is  ramp  scheduled  with  elevator  between  *) 
(*  -4.5  and  1.5  degrees  elevator  movement  and  a  constant  *> 

(*  +  or  -  0.2  outside  of  those  values.  Because  the  trim  *) 

(*  elevator  position  for  this  flight  condition  (350  KIAS,  *) 
(*  15,000  feet  PA)  is  approximately  -4.5  degrees,  the  ramp  *> 
(*  start  and  stop  points  have  been  shifted  to  reflect  a  trim  *> 
(*  elevator  position  at  -4.5  degrees.  *> 


IF  XC 193  <0.0  THEN 
k 1 : =1 . 0 
ELSE 

IF  X C 193  >  0.10471  "rHEN 
k 1 : =-l . 0 
ELSE 

k 1 : =-0. 5  -  (XC 193—0. 07853) $  19. 1 ; 

FC253 : =-XC253  +  0.25*Rdot; 
dr2:=(-0.2*kl*XC243) ; 

IF  ABS (drpi lot) <0.02  THEN 
Begin 

FC263 : =— 2*XC263  +  2*Ay; 

FC273 : =0. 0009*XC263 ; 

dr  1 : =0 . 003*  X  C  26  3  +  XC273; 

IF  drl  >  0.07  THEN  drl: =0.07 

ELSE  IF  drl  <  -0.07  THEN  drl: =-0.07; 

End; 

drcommand :=XC 25 3  *■  dr  1  +  dr2  drpilot; 

IF  drcommand >0. 24  THEN  drcommand : =0. 24 
ELSE 

IF  drcommand<-0. 24  THEN  drcommand : =— 0. 24; 
FC283  :  =-33.  3*XC283  33.  3*drcommand ; 

END; 


PROCEDURE  XPRIME (VAR  X:  VECTOR:  VAR  F:  VECTOR); 


BEGIN 

COEFFICIENTS  <  X , adot , Qbar , A1 pha, Beta , Cx , Cl , Cm, Cn , An , Ay) 
A7DYN (X, Qbar , Alpha, Cx , Cl , Cm, Cn , An , Ay, F, adot ) ; 

ASSIGN (F, Pdot , Qdot , Rdot ) ; 

TGTDYN (X, F, Lazimuth, Lei evati on , Antgt ) ; 

ELEVATOR ( X , An , Qdot , P I LOTel ev , F > ; 

AILERON ( X, PILOTai 1  ,F> ; 

RUDDER (X, PILOTrud, Ay, Rdot, F) ; 

END; 


PROCEDURE  PR INTOUT (VAR  T : REAL;  X  :  VECTOR); 

VAR  I  :  INTEGER; 

J  :  REAL; 
b  :  CHAR; 

F  :  INTERACTIVE: 

BEGIN 

WRITELN(T) ; 

J:=COUNT*Tint  -  0.01; 

IF  T< J  THEN  EXIT (PRINTOUT) 

ELSE 

IF  C0UNT=1  THEN 
BEGIN 

WITH  OUTPUT^  DO 

FOR  la  — 1  TO  30  DO  STATEC 13 : =0. 0; 

END; 

COUNT: =C0UNT+1; 

WITH  OUTPUT^  DO 
Begin 

FOR  I : =1  TO  30  DO  STATEC I  3 : =X C I  3 ; 
Normal accel : =An ; 

Tgtaccel :=Antgt; 
az i muth : =Laz i muth ; 
elevation: =Lel evati on; 

Sideaccel :=Ay; 

AOA: =A1 pha; 

Sidesl ip: =Beta; 

T i me: =T ; 

End; 

PUT (OUTPUT) ; 

END; 


PROCEDURE  RKDF ( VAR  XO: VECTOR;  TO, Tf i nal , Tdel ta : REAL  : 
VAR  Xfinal:  VECTOR); 


N :  I  N  t  p  ( 


<*  This  procedure  is  a  fourth-order  Runge-kutta 
(*  solver.  (Ref:  Numerical  Methods,  Robert  Hornb»i  >  . 
(*  It  requires  a  TYPE  declaration  in  the  main  <■< 
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(*  VECTOR  =  ARRAYC 1 . . 103  OF  REAL  *> 
<*  and  a  VAR  declaration  for  the  parameters  in  the  above  parameter  *) 
(*  list.  It  also  requires  a  procedure  XPRIME  and  PRINTOUT.  The  t) 
(*  XPRIME  procedure  contains  the  equations  to  be  solved.  The  t> 
(*  PRINTOUT  procedure  contains  the  output  format.  t) 


VAR  I,K,NN  :  INTEGER; 

X1,X2,X3,F0,F1,F2,F3:  VECTOR; 

Iterations,  Tend,  Tbidelta,  Thexdelta  :  REAL; 


BEGIN  (*  RKDF  *) 

FOR  I  s*  1  TO  30  DO 
F0CI3s*0.0; 

Tend  s=  TO; 

Tbidelta  :*  Tdelta/2; 

Thexdelta  :**  Tdelta/6; 

Iterations  s=  (Tfinal  -  T0)/Tdelta; 

NN  s«  ROUND ( Iterations) ; 

FOR  K  s-  1  TO  NN  DO 
BEGIN  (*  K  -  Loop  *) 

XPRIME (X0,F0); 

FOR  I  is  1  TO  N  DO 

X1CI3  s«*  X0CI3  +  (Tbidelta)  *F0C II; 
XPRIME (XI, FI ) ; 

FOR  I  s*  1  TO  N  DO 

X2C 1 3  s-  X0CI3  +  (Tbidelta) *F1CI3; 
XPRIME (X2,F2); 

FOR  I  :»  1  TO  N  DO 

X3CI3  :«  X0CI3  +  Tdel ta*F2C 13; 
XPRIME (X3,F3)  ; 

Tend  :=  Tend  +  Tdelta; 

FOR  I  :=  1  TO  N  DO 


Begin 

Xf inal C 1 3  X0CI3  +  (Thexdel ta) * (FOC 1 3  +•  2*(F1CI3  +  F2CI3)  + 

F3C 13); 

X0CI3  *  =  Xf inal C 1 3 ; 

End; 

PRINTOUT (Tend, Xf inal ) ; 

END;  (*  K  -  Loop  I) 

END;  (*  RKDF  *) 


BEGIN  (*  MAIN  PROGRAM 
(*  INITIALIZE  VARIABLES  *> 
FOR  II* =  1  TO  30  DO 
Begin 

FC 1 1  3  s  =0. 0; 

XOC 1 1 3  * “0.0; 

End; 
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■''**.  V, 


Qbars=0.0; 

Pdots “0.0; 

Qdots “0.0; 

Rdots“0.0; 

adot:“0.0; 

Antgt : =32. 2; 

Ptgts“0.0; 

PILOTelevs “0.0; 

PILOTai 1 : =0. 0; 

PIL0Truds=0.0; 

Lazimuth: “0.0; 

Lei  evati  on  s  “0. 0; 

Cx :  =0.0; 

Cl :  =0.0; 

Cm :  =0 . 0 ; 

Cn:=0.0; 

An: “0.0; 

Ays “0.0; 

Ns “28; 
bs='  '; 

COUNT: “1; 

(*  END  INITIALIZATION  *) 

WritelnC'  If  you  want  any  NONzero  state  values,  enter  them  one 
WritelnC'  at  a  time  when  prompted  by  first  entering  the  number 
WritelnC'  of  the  state,  a  space,  then  the  value.  ' ) ; 

Wr  i  tel  n  C '  Ceg.  X0C10J—4  would  be  10  -4<return>  '); 


WritelnC'  at  a  time  when  prompted  by  first  entering  the  number 
WritelnC'  of  the  state,  a  space,  then  the  value.  ' ) ; 

Wr  i  tel  n  C '  Ceg.  X0C10J—4  would  be  10  -4<return>  '); 

REPEAT 

WriteC'  Input  the  state  number  and  value  (zero  when  done)  ') 
Read (J) ; 

IF  <J<0>  OR  <J>N)  THEN  WritelnC'  Reenter  ’> 

ELSE 

IF  J»0  THEN  Writeln 
ELSE 

Readln (XOC JD) ; 

UNTIL  J“0; 

WriteC'  Enter  PILOTail  *  '); 

Readln CPILOTail); 

WriteC'  Enter  PILOTrud  “  '); 

Readln CPILOTrud) ; 

WriteC'  Enter  initial  time  =  '); 

Readln CTO) ; 

WriteC'  Enter  final  time  “  '); 

Readln  CTf inal ) ; 

WriteC'  Name  of  output  files  '); 

Readln  Cdi skf i 1 e) ; 

IF  LENGTH Cdi skf ile)“0  THEN  EXIT (PROGRAM) s 

IF  POSC' . ' ,diskfile)“0  THEN  di skf i 1 e: “CONCAT (di skf i le. '. DATA' > 
REWRITE (OUTPUT, di skf ile) ; 

RKDF(X0,T0,Tf inal, Tdelta,N,Xf inal > ; 

CLOSE (OUTPUT, LOCK) ; 

END.  (*  MAIN  PROGRAM  t) 


“  '); 


1 

• 

2*7  i 

v 

PROGRAM  PRINT; 

• 

i 

TYPE  VECTOR  « 

ARRAYC1. .303 

VAR  I, COUNT  : 

INTEGER; 

b  :  CHAR; 

disk-file  :  STRING; 

X  s  VECTOR; 

F  s  INTERACTIVE; 

INPUT  s  FILE  OF  RECORD 

STATE  :  VECTOR; 

Normal accel , Tgtaccel , azimuth, elevation, Si deaccel , 
AOA, Sideslip, Time  :  REAL; 

END; 


A 

1 


BEGIN 

Write ('Name  of  input  file:  *); 

Readln (diskf i le) ; 

IF  LENGTH (diskf ile)-0  THEN  EXIT (PROGRAM) ; 

IF  POS(*.*, diskf ile)-0  THEN 

diskf ile: -CONCAT (diskf ile, ' .DATA' ) ; 

COUNT: -1; 

RESET ( INPUT, di skf i 1 e) ; 

IF  COUNT-1  THEN 

BEGIN 

b:«'  '; 

REWRITE (F, 'PRINTER: ' ) ; 

Writeln(F,CHR(lS) ) ; 

Write (F, '  T  ',b:7, 'U',b:7, 'V',b:7,'W',b:7, 'P',b:7,'Q'>; 

Write (F,b:8, 'R' , b:6, 'PHI' ,b:3, ' THETA' ,b:3, 'PSI' ,b:6, 'de' ,b:6, 'da' 
Wri teln (F, 'dr' ,b:6, 'An' ,b:7, ' Ay' ,b:3, 'Beta' ,CHR(18) > ; 

CLOSE (F) ; 

END; 

REWRITE (F, 'PRINTER: ' ) ; 

Writeln(F,CHR(13) ) ; 

REPEAT 

IF  COUNT >1  THEN  GET (INPUT); 

WITH  INPUT'*  DO 
Begin 

FOR  I:-  1  TO  9  DO  XCI3:«STATECI3; 

XC 193 : — STATEC 193 ; 

XC24J:— STATEC243; 

XC283 : -STATEC283 ; 

T:-Time; 

An: -Normal accel ; 

Ay: -Si deaccel ; 

Beta: -Sideslip; 

End; 

COUNT: -COUNT* 1; 

Write(F,T:3:  1,  XC  1 3:8:  4,  XC23:8:5,  XC33:8:3,  XC43 : 8:5,  XC53  :8:S,  XC63 : 8: 
Write  (F,XC73:8:5,XC83:8:3,  XC93 : 8: 5,  XC  193 :  8:5,  XC243 : 8:3,  XC283 : 8:  5) ; 


,b:6>  ;| 


V  >  v  v  ",*  *„■»  v 


v. 
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Writeln  <F,  An:  8s  5,  Ay:  8: 5,  Beta:  8:  S) 
UNTIL  EOF (INPUT) ; 

Writeln (F,CHR(18) ) ; 

CLOSE (F) ; 

END. 


APPENDIX  C 


Development  of  Target  Dynamics  and  Relative  Kinematics 


DEVELOPMENT  OF  TARGET  DYNAMICS  AND  RELATIVE  KINEMATIC  EQUATIONS 

The  target  was  modelled  as  a  point  mass  with  a  velocity  and 
angular  orientations.  The  angular  orientations  were  necessary 
since  they  provide  important  feedback  cues  to  the  pilot  for  the 
air-to-air  tracking  task.  For  example,  the  pilot  needs  to  know 
the  bank  angle  of  the  target  so  that  he  can  match  his  bank  angle 
to  the  target  to  be  able  to  stay  in  the  target  turning  plane,  or 
at  least  in  a  parallel  turning  plane.  Actually,  for  the  tracking 
task,  the  pilot  only  needs  to  know  the  relative  difference  in 
angular  orientations  between  his  aircraft  and  the  target.  The 
absolute  orientation  is  of  use  for  air-to-air  tracking  "situation 
awareness"  (ie.  present  and  future  energy  states). 

The  following  equation  developments  for  the  target  dynamics 
and  relative  kinematics  are  provided  as  support  for  the 
discussions  in  the  body  of  this  report. 

TARGET  DYNAMICS 

The  first  step  in  describing  the  relative  motion  between  two 
bodies  in  motion,  is  to  reference  their  motion  to  a  fixed 
(inertial)  reference  frame.  For  the  air-to-air  tracking  task, 
the  earth  surface  fixed  reference  frame  will  suffice  as  an 
equivalent  inertial  frame.  The  target  can  be  referenced  to  the 
earth  fixed  frame  through  a  sequence  of  three  rotations. 


.<  v  "J»  ■> ■>  ■.»  •>»>  v -ji  v  v  v  v  i 


Assume 


el'e2*e3  =  orthogonal  earth  fixed  unit  vectors 
^1*^2* ^3  =  ort^°9onal  target  fixed  unit  vectors 

/vm 


V^o. 


A 

+  v. 


■>  e. 


A 

*e3 


VT  =  Vvl 


v3  is  normal  to  the  target  wings. 


The  three  rotations  will  be: 

A  /V  ^ 

(fj  about  e3  to  define  f^.f^f^^ 
9  about  ?2  to  define  g^g^^ 
6  about  ^  to  define  vi»v2/V3 


Top  View 


Side  View 


To  find  ATe: 

V  =  ed_  (VT)  =  ed_  (Vvx)  =  vd_  (VCX)  +  wve  X  (vOx) 
dt  dt  dt 

=  Vvx  +  (£vx  +  0g2  +  ^f3)  X  (Vvx) 

Using  the  relationships  defined  above  to  transform  the  g2  and  f 
terms  into  v  coordinates,  we  gets 

A^,e  =  Vv^  +  ($5v^  +  0cos0v2  -  ©sin^v^  -f'sinOv^  +  ^cosOsin^v 
+  J^cosOcos^Vj)  X  Vv^ 

*  V\?^  +  (V{fcos0cos^  -  V0sin^S)v2  +  (-V0cos^ 

-  VfcosOsin^) 


[Cev]  =  [Cef][Cf9][C9v] 


cos^cos©  -sinjt'cos©  +  sinGcosy/sinflS  sinf6sin/  + 

sin©cos/cos$£ 

cosGsin^  cos^Scos^  +  sin©siny/sinf6  -cos^sinj6  + 

sinGsin^cos^ 


-sin© 


cosOsin^ 


cos0cosf6 


Now  using  F  =  mA  to  get  the  target  equations  of  motion. 


mge^  =  mg(-sin0v^  +  cosGsin^^  +  cosGcos^v^ ) 

A  • 

:  mV  =  T  -  D  -  mgsin© 

Vj :  (V^cosGcos^  -  VGsin^)  =  gcosGsin^ 

$2 :  ~(V0cosf6  -  VfcosGsin^S)  =  gcosGcos^  -  F^/m 

assume  V  =  0,  $  =  Pfc  t  =  given,  Fa/m  =  given  (target  "g") 

•  A  • 

To  solve  for  ,  solve  the  Vj  equation  for  V©  and  substitute  into 
the  v^  equation.  Then  through  algebraic  manipulations: 
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tjj  =  F^/m  sin# 

V  cos9 

Substituting  this  back  into  the  V0  equation  and  solving  for  0 

Q  -  Fa/m  cos#  -  qcosQ 
V 

Relative  Kinematics 


Refering  to  figure  4,  the  vector  of  interest  for  the 
relative  kinematics  is  ?T/A  (the  position  of  the  target  with 
respect  to  the  attacker).  In  the  earth  fixed  frame: 

?T/A  "  xSl  +  yg2  +  z«3 

e±.  <?T/A>  -  V  -  V  -  Vl  -  (u£l  +  vfi2  +  Wb3 ) 

dt 


• 

X 

# 

y 

m 

[Cev] 

'VT 

0 

-  cceb] 

p  1 

u 

V 

e 

z 

e 

_0_ 

V 

w 

The  desired  quantity  is  the  motion  of  the  target  with 


respect  to  the  attacker.  Therefore,  describing  rT/A 
attacker  body  frame: 


in  the 


?T/A  =  rlbl  +  r2b2  +  r3b3 


(rm/J  =  bd  (r_ )  +  wbe  X  r„ 


—  '  T/A'  —  'tT/Ai 


Solving  for  bd  (*r/A): 


*4-  (fbT/A> 


^  <feT/AJ 


wbe  X  rb. 


*S_  '^T/A' 


V  e  -  V  e 
T  VA 


-be  -b 
w  X  r  „ 


Therefore: 


0  -R  Q  r 


R  0  -P  r 


-Q  P  O  r 


rx  0  -R  Q  r±  VT  U 

r2  -  -  R  0  -P  r2  +  [Cbe][Cev]  0  -  V 

r„  .  -Q  P  Or..  0  W 


The  direction  cosine  matrix  [Cev3  was  developed  earlier.  The 
matrix  [C  ]  can  be  developed  in  the  same  manner.  Using  the  same 
standard  sequence  of  rotations,  yaw,  pitch,  and  roll,  we  obtain 
three  direction  cosine  matrices  which  are  combined  to  arrive  at 
the  overall  direction  cosine  matrix  ([C  ]). 


[c  ]  = 


cosGcosf' 


sin^cosG 


-sinG 


cosf'sin^sinG  sin^sin^sinG  +  cos^cos(6  cosGsin>6 

-  sin^cos^ 

cos^cos^sinG  sin^/cos^sinG  -  sin^Scosp'  cosGcosj6 

+  sinj/sin^ 


[cbv]  =  [cbe][cev] 


C11  c12  c13 

=  C21  C22  C23 

C31  C32  C33 


It's  evident  from  the  vector  multiplier  to  [C  ],  that  only  the 
first  column  of  [Cbv]  is  of  interest.  The  subscript  T  is  used 
here  to  distinguish  target  angles  from  attacker  angles. 


C11  =  cos^c^pCosGrpCOsGcos tj/  +  cosG.pSin^sin^'cosG  +  sinGTsin  G 

C21  *  cos^TcosGTcos0/sin^sinG  -  cosf^cosGfpSin^cosjiS 

+  cosGTsin^Tsin^sin^sinG  +  cosG.j.sin^pCOsp'cos^ 

-  sin  6-cosGsin^ 


=  cos^Tcos©Tcos^sine  +  cos^Tcos©Tsiny'sin0 

+  coseTsin^Tsin(»,cos^sin0  -  sin/Tcos©Tcos^sinj6 
-  sinQ^cosQcosflJ 

Therefore: 


II 

H 

*U 

Rr2  "  Qr3 

+ 

C11VT 

-  U 

r2  = 

"Rrl  +  Pr3 

+ 

C21VT 

-  V 

CO 

II 

Qrl  -  Pr2 

+ 

C31VT 

-  w 

To  describe  how  the  target  would  appear  in  the  attacker  HUD 
"X_  =  azimuth  angle  of  target  in  HUD 

a 

■  arctan  (r2/r1) 

>e  ■  elevation  angle  of  target  in  HUD 


=  arctan  (-r-/r,  ) 


DEVELOPMENT  OF  LINEAR  EQUATIONS 


In  this  section  the  details  of  developing  the  linear  model 
of  the  air-to-air  lateral-directional  tracking  task  are  shown 
along  with  the  linear  mechanical  and  CAS  roll  control  system. 

The  starting  point  for  the  linear  model  of  the 
lateral-directional  air-to-air  tracking  task  is  the  linearized 
lateral-directional  equations  of  motion  for  the  aircraft.  The 
linearized  equations  were  obtained  from  Roskam  and  are  shown 
below. 

*  ( l/U )  A  +  g *5  +  Ypp)  +  (Yr/U  -  l)r 

+  (l/U)(YjaXa  +  YSrSr) 

P  **  LJj  /*  +  L'pp  +  L*rr  +  L*£a£a  +  L*jrJr 
r  =»  N^  A  +  N'pP  +  N'rr  +  N 1  £a  Ja  +  N'Jrfr 
b  -  p 


where 

L' 


C(La  +  <Ixz/Ixx,V  1  (l  -  <(Ixz> 


N'  -  [ (N  +  (I  /I  )L  )  /  (1  -  ((I  )  /I  I  ))] 

a  ‘•'a  xz'  zz  a  '  xz  xx  zz 


Arranging  these  equations  in  state  variable  form  (x  =  [A]x  + 
[B]u)  where  xT  *  t/3  p  r  56]  and  uT  *  [£_  S_]  we  obtain  the  A 

Q  1 


matrix: 


and 


Yia/° 

Y&/U 

L’<Ta 

L'Sr 

NVa 

N7r 

0 

0 

The  stability  parameters  were  obtained  from  an  LTV  report 
(Reference  6)  and  are  shown  in  Table  Dl.  Substituting  in  these 
values  we  get: 


[A]  =  I  - 


0.162 

26.23 

4.546 

0 


0.00089 

-3.008 

0.05664 

1 


-0.9979 

0.9596 

-0.5298 

0 


0.05079 

0 

0 

0 


0.01132 

34.55 

0.0614 

0 


0.0453 

5.972 

-5.307 

0 


Changing  the  state  variable  /3  to  v  where  v  */^U,  we  get  a  new 
state  vector,  x  ■  [v  p  r  and  new  plant  and  control  matrices 
which  differ  from  the  previous  matrices  in  that  the  first  row  of 
both  matrices  is  multiplied  by  U  and  the  first  column  of  the 
plant  matrix  is  divided  by  U.  This  was  done  primarily  as  a 
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TABLE  D1 


Lateral-Directional  Stability  Parameters 


Yft 

(ft/sec2) 

-102.6961 

YP 

( ft/sec) 

0.5643 

Yr 

(ft/sec) 

1.3516 

YSa 

(ft-rad/sec2) 

-7.1836 

YSr 

(£t-rad/sec2) 

28.7641 

(sec-2) 

-25.7350 

L 

(sec-1) 

-3.00146 

P 

Lr 

(sec-1) 

0.9022 

(rad/sec) 

34.3848 

(rad/sec) 

5.4094 

(sec-2) 

3.9938 

N 

P 

(sec-1) 

-0.0067 

Nr 

(sec-1) 

-0.5096 

Nfa 

(rad/sec) 

0.7854 

"jr 

(rad/sec) 

-5.0721 

cosmetic  change  because  the  relative  kinematic  equation  developed 
earlier  contained  a  v  term. 

Next  assume  no  rudder  inputs  to  eliminate  and  add  a  roll 
attitude  pilot  model.  The  effect  of  this  roll  attitude  pilot 
model  is  to  couple  the  roll  and  spiral  poles  together  and  make 
them  a  complex  conjugate  pair  with  good  damping  and  a  fast 
response  time.  The  roll  attitude  pilot  model  is  realized  by 
first  closing  an  inner  roll  rate  feedback  loop  and  then  using 
roll  attitude  feedback  as  the  outer  loop.  This  is  illustrated  in 
figure  Dl.  Another  effect  of  using  roll  attitude  feedback  is 
changing  the  control  variable  from  f  to  6. 

In  figure  Dl,  the  transfer  function  p/$  is  obtained  from 
the  A  and  B  matrices  shown  earlier. 

£  -  34.55  s  (s  +  0.351  +  j2.13) _ 

£  (s  +  0.036) (s  +  0.34  +  j  2 . 1 ) (s  +  3) 

a  •“ 

To  simplify,  assume  that  the  spiral  effect  is  eliminated  by  the 
zero  at  the  origin  and  that  the  dutch  roll  is  not  present  due  to 
the  zeroes  close  by.  This  implies  that  the  spiral  and  dutch  roll 
are  not  prevalent  in  an  aileron  roll  and  that  the  response  is 
primarily  due  to  the  roll  mode  root.  This  is  true  for  most 
aircraft  and  is  also  true  in  the  case  of  the  A-7D.  Thus,  with 
the  above  simplifications: 

£  -  35,  _ 

£  s  +  3 
a 

The  values  for  Kp  and  were  chosen  to  give  good  response 
and  adequate  damping.  By  choosing  Kp  *  0.1  and  =  0.5,  we  get 
the  closed  loop  transfer  function: 


s  +  6.5s  +  17 . 5 


This  yields  a  final  value  of  one,  a  natural  frequency  of 
approximately  four  radians/second,  and  a  damping  of  approximately 
0.77  which  are  within  the  average  pilot's  abilities. 

The  inclusion  of  this  roll  attitude  pilot  model  into  the 
linear  model  was  done  by  first  writing  the  old  control  variable 
(S _)  in  terms  of  the  new  control  variable  (d  ). 

a  C 

*  0.5^c  -  0.5*S  -  O.lp 

Next  the  feedback  matrix  (K)  is  used  to  obtain  a  new  plant 
matrix. 

[K]  «  [0  0.1  0  0.5] 

[A]  -  [A]  -  [B][K] 


-0.162 


1.282 


-632.7 


35.79 


-0.0414  -6.463 


0.9596  -17.28 


0.0072 


0.0505  -0.5298  -0.0307 


Note  that  here  xT  «  [v  p  r  ^]. 
u*  ■  4>c,  therefore, 

CB']  -  0 . 5[B] 


[B‘] 


-3.588 


17.28 


0.0307 


The  next  step  was  the  addition  of  a  yaw  SAS  to  increase  the 
dutch  roll  damping  (ie.  move  the  dutch  roll  poles  away  from  the 


SB* 


m 


imaginary  axis).  To  accomplish  this,  the  yaw  damping  term  (Nr‘, 
the  3,3  element  of  [A])  was  artificially  increased  to  -2.0.  The 
effect  was  to  increase  the  dutch  roll  damping  from  0.16  to  0.50. 

The  modified  linear  aircraft  model  was  then  augmented  with  a 
heading  term  (^) .  The  heading  term  was  included  since 
lateral-directional  tracking  involves  lateral  pointing  which  can 
be  represented  with  heading  angle. 
iy  -  r 

To  this  fifth  order  linear  system  was  added  a  lateral 
relative  kinematic  term  which  represents  the  lateral  tracking 
error.  The  equation  used  was  a  linearized  version  of  the  lateral 
element  of  the  relative  kinematic  equations. 

r2  -  -rrx  +  pr3  +  VT(-^)  -  v 

Me  now  have  a  sixth  order  state  variable  representation  of 
the  lateral-directional  air-to-air  tracking  task. 


£  v  p  r  ^  r,  ] 


u  -  4>r 


-0.162  1.282 


-632.7  35.79 


-0.0414  -6.463  0.9596  -17.28 


0.0072  0.0505 


-0.0307 


-37.5 


-1500 


Note:  CA]  is  for  147  mil  roll  axis  depression. 


CB] 


-3.585 

17.28 

0.0307 


0 


0 


0 


The  linear  model  of  the  air-to-air  tracking  task  does  not 
contain  any  of  the  control  system  elements  which,  from  the 
nonlinear  analysis,  have  an  effect  on  aircraft  response  to  pilot 
inputs.  The  control  systems  could  not  be  incorporated  into  the 
above  linear  model  because  of  limitations  on  the  size  of  state 
variable  models  in  TOTAL.  In  order  to  compare  the  linear  fourth 
order  lateral-directional  aircraft  model  to  flight  test,  the 
linearized  mechanical  and  CAS  roll  axis  controls  were  added  to 
the  aircraft  model.  Referencing  figure  2,  the  equations  for  the 
mechanical  and  CAS  roll  axis  control  are: 

A.  «  -3.33A.  +  0.7F 
1  1  sa 

A2  «  10AX  -  10A2 

A,  -  -12.8A-  +  2F„, 

J  3  sa 

A4  -  1.238A3  -  12.5A4 
fac  ■  o.i a2  -  O.lp  +  a4 
fa  -  2A2  -  2p  +  20A4  -  20  £a 
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The  primary  nonlinear  element  that  was  neglected  in  this 
representation  was  the  series  servo  limit  on  the  CAS  authority. 
The  new  state  vector  is  xT  =  [  v  p  r  ^  £ a  3 ,  and  the 

new  control  variable  is  lateral  stick  force  (Fo:i).  The  A  matrix 
is  shown  on  the  next  page.  The  B  matrix  is: 

[B]T  =  [0000  0.7  02003 

This  is  the  ninth  order  linearized  roll  CAS  model  that  was 
compared  to  flight  test  for  p/F  frequency  response. 
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